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Abstract 


This  project  aims  to  provide  some  insight  on  a  local  company’s  (s,  S)  inventory 
control  system.  It  does  so  by  providing  three  completely  new  ( s ,  S)  approaches  and 
comparing  their  simulated  performance  to  the  current  system.  Recommendations  are 
made  based  on  the  results. 


1  Introduction:  Quantegy,  Inc. 

Quantegy,  Inc.  (formerly  known  as  Ampex)  is  an  Opelika  company  that  manufactures  mag¬ 
netic  tape  for  worldwide  distribution.  It  has  been  the  sole  producer  of  magnetic  tape  in  the 
United  States  since  last  year,  when  Quantegy’s  only  competitor,  3M,  left  the  market.  The 
company  has  been  in  business  for  over  50  years.  Annual  sales  presently  average  $120  million. 
Management  estimates  that  nearly  $15  million  is  currently  invested  in  worldwide  inventory. 

The  facility  in  Opelika  is  Quantegy’s  only  manufacturing  plant.  It  also  serves  as  the 
company’s  headquarters.  Quantegy  maintains  ten  regional  distribution  centers  in  North 
America,  six  in  Europe,  one  in  Asia,  and  one  in  Australia.  These  regional  facilities  conduct 
Quantegy’s  business  in  their  respective  areas,  having  their  inventories  periodically  replen¬ 
ished  with  shipments  from  Opelika. 

Quantegy  uses  an  (s,  S )  inventory  control  policy  to  keep  its  regional  facilities  stocked. 
Specifically  the  inventory  at  each  facility  is  reviewed  from  one  to  three  times  per  week. 
When  the  inventory  level  for  a  particular  stock  keeping  unit  (hereafter  an  SKU)  is  found  to 
be  below  its  reorder  point  s,  then  a  replenishment  order  is  placed  to  Opelika.  The  amount 
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of  the  order  is  equal  to  the  difference  between  the  current  inventory  level  and  that  SKU’s 
order-up-to  level  S.  All  order  amounts  are  integer  multiples  of  the  standard  package  size  for 
that  particular  SKU.  The  distribution  center  can  expect  to  receive  the  order  after  a  delivery 
lagtime. 

The  reorder  point  s  and  order-up-to  level  S  change  from  week  to  week,  based  on  the 
demand  forecast  for  that  particular  SKU.  The  forecast  for  the  next  week’s  demand  is  simply 
the  average  demand  over  the  last  13  weeks  (i.e.  a  13  week  moving  average).  The  order-up-to 
level  S  is  this  forecast  multiplied  by  a  number  called  the  ABC  factor  and  then  rounded 
to  the  nearest  standard  package  size.  Quantegy  classifies  its  SKU’s  according  to  their  own 
ABC  system,  where  A  items  are  the  most  important.  Their  system  does  not,  however,  use 
the  classic  Pareto  rule  (80/20)  [1].  Instead  each  sales  region  requests  a  specific  classification 
that  is  subject  to  the  main  planning  department’s  final  decision.  The  ABC  factor  is  just  a 
number  that  management  creates  and  can  adjust  in  times  of  irregular  demand.  The  factor 
is  typically  larger  for  A  items  (e.g.  3)  than  for  B  or  C  items  (e.g.  2.5).  The  reorder  point  s 
is  90%  of  the  order-up-to  level  S,  again  rounded  to  the  nearest  standard  package  size. 

For  example,  assume  a  B  item.  Its  ABC  factor  is,  arbitrarily,  2.5.  Assume  the  average 
demand  for  the  product  over  the  last  13  weeks  is  26.78  SKU’s.  The  standard  package  size 
for  this  product  is  10.  It  follows  then  that  the  order-up-to  level  S  is  26.78*2.5  =  66.95  «  70. 
The  reorder  point  s  is  90%  *  66.95  =  60.255  «  60. 

Each  region  also  provides  Opelika  with  its  own  demand  forecasts  to  take  into  consider¬ 
ation  when  evaluating  its  inventory.  The  general  idea  is  that  the  regional  personnel  should 
know  when  something  irregular  is  happening  in  their  area  of  responsibility.  But  these  fore¬ 
casts,  however,  are  notoriously  large.  Regional  managers  evidently  provide  inflated  numbers 
to  maximize  their  own  inventories  at  the  expense  of  the  other  regions.  As  a  result,  the 
regional  forecasts  are  virtually  ignored  by  Opelika  except  in  the  most  extreme  cases. 

Quantegy  measures  its  fillrates  by  calculating  the  percentage  of  orders  that  are  fulfilled 
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by  on-hand  demand.  The  sizes  of  the  orders  and  the  stockouts  are  ignored.  Management 
typically  expects  A  items  to  achieve  a  95%  fillrate  and  B  items  to  achieve  a  90%  fillrate.  C 
and  D  items  are  not  expected  to  achieve  anything  close  to  90%. 

2  Project  Overview 

The  project  objective  is  to  evaluate  Quantegy’s  current  (s,  S )  inventory  control  system 
(section  6)  against  the  three  different  (s,  S )  approaches  presented  in  this  paper  (sections 
7-9).  Evaluation  comes  in  the  form  of  a  deterministic  simulation  (section  5)  that  uses  actual 
demand  data  provided  by  analysts  in  Opelika  (section  4).  The  results  are  then  compared  to 
the  present  system’s  simulated  performance  and  recommendations  are  made  (section  10). 

The  methodology  used  is  relatively  straightforward.  The  data  is  first  partitioned  into 
two  sets.  The  first  set,  called  the  sample,  is  used  to  initialize  the  approach  being  looked 
at.  (A  “policy”  or  “approach”  is  simply  how  the  s  and  S  values  are  being  computed.) 
No  modifications  to  the  control  policy  are  allowed  after  initialization.  The  policy  is  then 
evaluated  via  simulation  on  the  second  set  of  demands,  called  the  ex  post  [2].  Assuming  a 
valid  simulation,  this  approach  discovers  what  would  have  happened  (i.e.  fillrates,  number 
of  orders,  average  on-hand  inventory)  had  Quantegy  actually  implemented  the  policy  in 
question,  during  the  times  represented  by  the  ex  post  set,  without  any  knowledge  of  the 
future  (ex  post)  demands. 

This  analysis  will  hopefully  provide  management  at  Quantegy  with  a  different  perspective 
and  some  fresh  insight  on  their  inventory  control  system.  This  project  does  not  claim  to 
offer  a  panacea  for  all  of  Quantegy’s  inventory  problems.  Instead  it  offers  some  practical 
alternatives  that  appear  to  have  significant  and  realistic  improvements  on  paper. 
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3  Literature  Review 


The  main  research  paper  behind  this  project  is  [14].  It  is  an  efficient  approach  that  uses 
renewal  analysis  to  find  the  minimum  cost  (s,  S )  policy  (long-run  expected  holding,  setup, 
and  penalty).  The  authors  themselves  improve  upon  the  work  of  [12],  who  solve  the  same 
problem  by  inefficiently  enumerating  all  solutions  in  the  (s,  A)  plane  and  selecting  the  policy 
with  the  minimum  cost. 

Another  paper  worth  mentioning  is  [7],  which  details  a  non-exact  heuristic  to  solve  the 
(s,  S)  problem  subject  to  a  Type  1  fillrate  constraint  (fraction  of  non-stockout  time  periods). 
This  project  attempts  to  do  the  same,  but  with  renewal  analysis  and  a  Type  2  fillrate 
constraint  (fraction  of  non-stockout  orders).  A  classic  paper  on  the  (s,  S)  problem  that  was 
looked  at  but  not  used  is  [9],  where  the  author  introduces  three  methods  for  computing 
approximately  optimal  (s,  S)  policies  for  discrete  time  shortage  cost  systems. 

[13]  provides  Croston’s  method,  an  advanced  forecasting  technique  used  in  the  project. 
[4]  enables  the  avoidance  of  the  common  mistake  of  underestimating  fillrates.  Numerous 
textbooks  provide  miscellaneous  support. 

4  Data 

Quantegy  was  able  to  provide  data  on  14  A  items  and  11  B  items  in  their  New  York  distribu¬ 
tion  center.  Specifically  they  provided  each  item’s  standard  cost,  and  the  quantity  shipped 
of  that  item  to  New  York  from  Opelika  every  week,  from  June  1993  to  October  1996  (172 
weeks  total).  Needless  to  say,  this  data  is  not  the  actual  demand  for  those  items.  It  is, 
however,  a  reasonably  close  substitute  and  the  only  information  Quantegy  could  give. 

The  raw  data  (attachment  1)  can  be  generally  described  as  erratic.  Each  item  displays 
its  own  behavior  such  as  trends,  bulges,  cutoffs,  and  long  periods  of  inactivity.  Graphs  of 
the  raw  weekly  demand  for  items  A14  and  B04  are  included  in  attachment  2  as  arbitrary 
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examples.  The  most  notable  characteristic  common  to  all  is  the  existence  of  outliers.  Each 
item  had  at  least  one  week  in  which  the  recorded  demand  was  more  than  three  standard 
deviations  from  the  mean  demand  [6].  Such  behavior  is  sure  to  pollute  even  the  most  robust 
control  policy’s  performance.  It  is  reasonable  to  assume  that  management  knows  about  such 
aberrant  orders  beforehand  and  can  prepare  accordingly.  To  reflect  this  assumption  and 
avoid  the  degradation  of  simulated  performance,  each  data  point  that  lies  more  than  three 
standard  deviations  from  the  mean  is  replaced  with  the  original  sample  mean.  A  summary 
of  the  exact  changes  is  included  in  attachment  3.  This  outlier  adjusted  data  (attachment  4) 
is  much  more  tractible  and  predictably  has  lower  means  and  standard  deviations. 

The  next  major  modification  dealt  with  changing  the  frequency  of  observations  from 
every  week  to  twice  per  week.  Such  a  change  is  necessary  because  the  inventory  levels  are 
reviewed  twice  per  week  at  the  New  York  facility.  Converting,  the  weekly  numbers  involves 
randomly  splitting  each  observation  into  two  half-week  observations.  But  some  extra  analysis 
(below)  is  needed  to  avoid  the  high  probablilty  of  these  half-week  observations  always  coming 
in  pairs,  an  outcome  that  was  pointed  out  at  a  presentation  of  this  project  at  a  local  APICS 
meeting  in  February  1997. 

Let  /  be  the  probability  of  demand  in  a  particular  week.  Estimate  it  by  dividing  the 
number  of  non-zero  observations  by  the  total  number  of  observations  (172).  Let  p  be  the 
probability  of  demand  in  a  particular  half-week  and  assume  that  half-week  demand  is  inde¬ 
pendent  and  identically  distributed.  The  probability  then  of  not  incurring  a  demand  in  any 
week  (1  —  /)  should  be  equal  to  the  probability  of  not  incurring  a  demand  in  both  the  first 
half  and  second  half  of  any  week  ((1  —  p)  *  (1  —  p)).  Solving  for  p  results  in  the  following 
equation. 


p  =  l-y/l-f  (1) 

[The  above  equation  can  be  generalized  for  any  number  n  of  time  periods  within  the  week 
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with  the  expression  p  =  1  —  l/l  —  /.] 

There  are  two  possible  outcomes  regarding  the  half-week  demands,  given  that  a  non-zero 
demand  has  occurred.  Either  (1)  the  demand  is  split  between  both  half-weeks  or  (2)  the 
entire  amount  occurred  in  either  the  first  or  second  half  of  the  week.  The  probabilities  of 
each  of  these  occurrences  can  be  calculated  using  conditional  probabilities  and  the  value  of 
p  above. 

•  P  (Demand  in  both  halves  |  Non-zero  weekly  demand) = p2 +2^  1  ^ 

•  P(Demand  in  only  one  half  |  Non-zero  weekly  demand) = 

The  conversion  procedure  logically  follows,  where  D  is  the  observed  weekly  demand,  Di 
is  the  resulting  demand  in  the  first  half  of  the  week  while  D2  is  the  resulting  demand  in  the 
second  half  of  the  week. 

Week  to  Half-week  Conversion  Procedure 

For  every  item 

•  Calculate  p 

•  For  every  week 

—  If  D  >  0  then 

*  Generate  x  =U(0,1) 

*  If  x  <  7?-  then  demand  in  both  halves 

•  Generate  y  =U(0,1) 

•  Di  =  y  *  D  rounded  to  the  nearest  standard  package  size 

•  D2  =  D  —  D\ 

*  Else  if  x  <  then  demand  in  first  half  only 

•  Di  =  D 

•  £>2  =  0 

*  Else  demand  in  second  half  only 

•  Dx  =  0 

•  D2  =  D 

—  Else  no  demand  in  either  half 
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*  D\  —  D2  —  0 


The  set  of  half-week  data  used  throughout  this  project  that  results  from  the  above 
stochastic  procedure  is  included  in  attachment  5. 

The  final  modification,  as  described  above,  is  to  partition  the  data  into  the  sample  and 
ex  post  sets.  The  last  year  of  data  (51  weeks  as  defined  by  Quantegy)  seemed  an  intuitively 
attractive  amount  of  time  for  the  ex  post  set,  which  provides  ample  data  for  initialization 
within  the  sample  (121  weeks).  This  partition  is  represented  by  the  horizontal  line  in  at¬ 
tachment  5. 


5  Simulation 


As  already  stated,  the  simulation  used  to  evaluate  the  inventory  control  policies  is  determin¬ 
istic.  It  simply  uses  the  half-week  demands  obtained  above  as  the  actual  demands  in  the 
simulation,  in  the  same  order  and  with  no  changes. 

Define  the  following.  .  . 
t  — >  the  time  period  in  half-weeks 

(e.g.  half-weeks  30.5  and  31.0  are  the  first  and  second  halves  of  week  31,  respectively) 

Yt  — >  the  demand  in  t 

st  — >•  the  reorder  point  in  t 

St  — »  the  order-up-to  level  in  t 

It  the  inventory  position  at  the  end  of  t 

Nt  — »  the  net  inventory  at  the  end  of  t 

Qt  — *  the  amount  of  the  order  placed  in  t 

qt  — »  whether  or  not  an  order  was  placed  in  t  (1  is  yes,  0  is  no) 

Ht  the  on-hand  inventory  at  the  end  of  t 
Bt  — >  the  amount  of  the  backorder  in  t 

The  simulation  then  is  just  the  act  of  making  the  following  calculations  for  every  t.  .  . 

It  =  It- 0.5  ~Yt  +  Qt- 0.5 
Nt  =  Nt- 0.5  —  Yt  +  Qt- 1.0 
Qt  =  St  -  It  if  It  <  st,  otherwise  Qt  =  0 
qt  =  1  if  Qt  >  0,  otherwise  qt  =  0 
Ht  =  Nt  if  Nt  >  0,  otherwise  Ht  =  0 
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Bt  =  —Nt  —  max( 0, 5t_0. 5  —  Qt~i.o)  if  Nt  <  0,  otherwise  Bt  =  0 

To  initialize,  Y0  =  0,  I0  =  N0  =  S0.5,  and  JV0.5  =  N0  -  Y0. 5. 

Some  things  are  worth  mentioning  here.  Quantegy  updates  the  s  and  S  parameters  every 
week.  As  a  result,  st  =  St+0.5  whenever  t  +  0.5  is  an  integer.  The  same  holds  true  for  S.  The 
equation  for  inventory  position  adds  in  the  order  from  the  last  period  because  the  delivery 
lag  to  New  York  is  half  a  week  long.  The  equation  for  net  inventory  adds  in  the  order  from 
an  entire  week  prior  to  keep  the  on-hand  inventory  as  accurate  as  possible  for  holding  cost 
estimates.  The  expression  for  backorders  prevents  previous  stockouts  from  being  double 
counted  if  a  demand  is  larger  than  normal. 

The  first  measure  of  performance  is  the  number  of  orders  made  by  a  policy.  The  fewer 
orders  a  policy  has  to  make,  the  better  in  terms  of  ordering  costs.  It  is  calculated  simply  by 
summing  the  qt  values.  The  next  measure  is  the  average  on-hand  inventory  and  is  calculated 
by  averaging  the  Ht  values  in  the  simulation.  Again,  this  measure  is  desired  to  be  as  small 
as  possible  to  cut  holding  costs.  The  final  measure  of  performance  is  the  fillrate,  calculated 
by  the  expression  below. 


EtBt 

EtY 


(2) 


This  is  commonly  referred  to  as  the  (3,  or  Type  2,  fillrate  and  is  different  from  the  fillrate 
that  Quantegy  uses,  commonly  called  the  a,  or  Type  1,  fillrate  [8].  The  difference  is  that 
the  /3  value  measures  the  units  of  demand  that  were  filled  with  on-hand  inventory  while  the 
a  value  measures  the  number  of  orders  (or  time  periods)  that  were  filled  entirely  by  on-hand 
inventory,  regardless  of  the  size  of  the  stockout.  The  j3  value  is  used  because  it  gives  better 
information.  Take  the  case  where  99  of  100  orders  were  filled  by  on-hand  inventory.  The  a 
fillrate  here  is  equal  to  99%  because  99  of  the  100  orders  were  filled.  But  now  assume  that 
there  were  10000  total  units  demanded  during  all  of  these  orders  and  that  the  one  missed  ■ 
order  (entirely!)  was  for  9901  units.  The  corresponding  (3  value  'here  is  1%,  a  much  more 
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indicative  number  for  the  situation.  Unlike  the  other  measures  of  performance,  the  fillrate  is 
desired  to  be  as  close  to  1  as  possible  (maximized)  and  no  smaller  than  the  service  objectives 
set  by  management. 


6  Quantegy’s  Current  Policy 

Recall  that  Quantegy  uses  a  dynamic  inventory  control  policy  where  the  values  of  s  and  S 

are  updated  every  week. 

Define  the  following.  .  . 

Yt  — >■  the  demand  forecast  for  t 

S°  -4  the  order-up-to  level  in  t  before  rounding 

s®  — »  the  reorder  point  in  t  before  rounding 

k  -4  the  ABC  multiplier  (3  for  A  items,  2.5  for  B  items) 

Then  for  every  integer  value  of  t  +  0.5.  .  . 

Yt  +  ut+0.5  =  Yt 

S°  =  k*  {Yt  +  lt+0.5) 
s°  =  0.9  *  S° 

St  =  St+0.5  —  rounded  to  the  nearest  standard  package  size 
st  =  St+0.5  —  s°  rounded  to  the  nearest  standard  package  size 

This  policy  uses  only  the  last  13  weeks  of  the  sample  set  to  initialize  the  moving  average. 
The  detailed  results  of  the  ex  post  simulations  for  all  25  items  are  included  in  attachment 
6.  The  table  below  summarizes  the  simulated  performance  of  Quantegy’s  current  inventory 
control  policy.  Overall  the  policy  made  a  total  of  1200  orders  for  all  25  items  (out  of  a 
possible  2550).  It  kept  an  estimated  $98,408.33  worth  of  inventory  on  the  shelves.  Of  the  25 
items,  15  had  fillrates  of  at  least  90%.  The  average  fillrate  for  the  entire  policy  was  86.8%. 

Since  Quantegy  does  not  keep  information  on  past  inventory  performance,  it  is  nearly 
impossible  to  exhaustively  validate  the  simulation.  Management  did,  however,  indicate  that 
the  numbers  were  “reasonable”  and  “seemed  right”  in  a  meeting  back  in  December  of  1996. 
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Performance  Measure 

Current  Policy 

Total  Orders 

1200 

Avg  OH  Inv  ($) 

98,408.33 

#  Fillrates  >  90% 

15 

Avg  Fillrate 

86.8% 

Table  1:  Performance  of  the  Current  Policy 

7  Alternative  1:  Different  Forecasts 

This  approach  tests  the  intuition  that  Quantegy’s  inventory  control  policy  would  perform 
better  with  a  more  accurate  and  sophisticated  forecast  than  a  13- week  moving  average.  It 
does  so  by  replacing  the  moving  average  with  three  alternatives.  No  other  changes  are  made 
to  the  policy. 

The  first  forecasting  alternative  is  the  naive  forecast,  where  the  demand  for  the  current 
week  is  simply  the  observation  in  the  previous  week.  It  uses  only  the  last  week  of  demand 
in  the  sample  set  to  initialize.  This  forecast  should  perform  poorly  compared  to  the  moving 
average  and,  as  a  result,  have  worse  inventory  performance.  It  is  defined  by  the  following 
equation,  where  t  +  0.5  is  again  an  integer. 

Yt  +  Yt+ 0.5  =  Yt- i.o  +  Pt— o.5  (3) 

The  second  alternative  is  to  use  the  mean  demand  over  the  entire  sample  set  as  the 
forecast  for  every  demand  in  the  ex  post  set.  Needless  to  say  it  needs  the  entire  sample  set, 
to  initialize.  This  forecast  results  in  a  static  policy  where  the  values  for  s  and  S  are  constant 
over  the  entire  simulation.  Like  above,  this  forecast  should  perform  poorly  both  in  terms 
of  accuracy  and  inventory  performance  compared  to  the  moving  average.  This  forecast  is 
defined  by  the  following  equation,  where  t  +  0.5  is  integer. 


Yt  +  Yt+ o.5  =  ^geyi  (4) 

The  final  alternative  is  called  Croston’s  method  [13].  It  is  a  variant  of  single  exponential 

smoothing  that  was  designed  to  forecast  random  and  intermittant  demand  like  the  data 

in  this  project.  Specifically  it  makes  separate  exponential  smoothing  estimates  for  both  the 

average  size  of  a  demand  and  the  average  interval  between  demands.  The  forecast  is  the  ratio 

of  these  two  values,  or  the  expected  demand  per  interval.  It  is  only  updated  when  a  demand 

occurs.  If  demand  is  occurring  in  every  time  period,  then  Croston’s  method  is  equivalent 

to  single  exponential  smoothing.  Unlike  the  other  alternatives,  this  forecast  should  be  more 

accurate  than  a  moving  average  and,  as  a  result,  have  superior  inventory  performance. 

Define  the  following.  .  . 

Zt  ->  the  estimate  for  the  size  of  demand  in  t 

Pt  — »  the  estimate  for  the  time  interval  between  demands  in  t 

d  — >■  an  integer  counter 

a  — >  smoothing  parameter  between  0  and  1 

Croston’s  method  follows,  where  t  +  0.5  is  an  integer. 

•  If  Yt  +  lt+o.5  =  0  then 


Initialization  here  uses  the  entire  sample  set  and  is  a  bit  more  involved.  Z- 0.5  +  Zq,  or 
the  initial  estimate  for  the  size  of  the  demand,  is  the  average  size  of  the  non-zero  demands 
from  y0.5  to  Tjo.o,  the  first  ten  weeks.  P-0.5  +  Po-,  or  the  initial  estimate  for  the  time  interval, 
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is  the  average  size  of  the  time  intervals  between  Y0. 5  to  F10.o,  cutting  off  the  last  interval  if 
undefined  by  time  10.0.  The  first  d  value  is  1.  The  initial  cm  value  is  0.1.  A  search  procedure 
(Excel  solver)  is  then  performed  on  a  to  find  a  local  minimum  for  the  mean  squared  error 
(hereafter  MSE)  of  this  forecast  over  the  sample  set.  Recall  that  the  MSE  is  simply  the 
average  squared  difference  between  forecast  and  observation  over  all  relevant  time  periods. 
The  initialization  results  are  included  in  attachment  7. 

Using  MSE  as  the  only  forecast  accuracy  measure,  the  initialized  naive  and  mean  forecasts 
are  less  accurate  over  the  ex  post  sample,  as  expected.  The  naive  forecast  results  in  MSE’s 
that  are  93%  more,  on  average,  than  the  MSE’s  produced  by  the  13  week  moving  average. 
The  mean  forecast  is  not  as  bad,  having  MSE’s  that  are  only  38%  greater.  Croston’s  method, 
on  the  other  hand,  is  slightly  more  accuarate  than  the  moving  average.  It  produces  MSE’s 
that  are  overall  1%  lower.  The  specific  accuracies  for  each  item  and  each  forecast  are  included 
in  attachment  8. 

The  simulated  performances  of  the  three  forecasts  are  summarized  in  the  table  below. 
The  complete  simulation  results  are  included  in  attachment  9. 


Performance  Measure 

Naive 

Mean 

Croston 

Total  Orders 

590 

1323 

1151 

Avg  OH  Inv  ($) 

173,144.69 

41,052.26 

98,448.15 

#  Fillrates  >  90% 

22 

6 

18 

Avg  Fillrate 

95.8% 

70.2% 

90.6% 

Table  2:  Performance  Comparison  For  Different  Forecasts 


8  Alternative  2:  A  Two- Stage  Heuristic 

This  approach  relies  heavily  on  the  renewal  analysis  presented  in  [14]  and  is  a  complete 
change  from  Quantegy’s  current  policy.  The  first  stage  involves  selecting  a  S  -  s  value  for 
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each  item  by  matching  its  classic  economic  order  quantity  (EOQ)  to  the  corresponding  order 
quantity  defined  by  the  renewal  equations.  The  second  stage  involves  increasing  each  item’s 
s  from  0,  and  likewise  S  from  S  -  s,  until  the  desired  fillrate  is  achieved  in  the  renewal 
analysis.  The  simulated  performance  of  these  static  s  and  S  values  is  then  evaluated  for 
each  item  in  the  ex  post  simulation  for  comparison  with  the  other  methods. 

Some  initialization  is  required,  however,  before  beginning  the  heuristic.  The  first  thing 
to  be  done  is  to  choose,  for  each  item,  a  discrete  probability  distribution  to  represent  its 


L./rh 


^^degaajxdj  This  was  accomplished  by  fitting  the  first  two  moments  of  the  half- week  demands 
in  the  sample  set,  divided  by  their  standard  package  sizes,  to  the  negative  binomial  distribu- 


tion  defined  below.  ^The  choice  of  distribution  is  appropriate  (as  compared  to  choosing  the 
traditional  normal,  poisson,  or  binomial  distributions)  since  the  coefficient  of  variation  (— ) 


for  each  item’s  half- week  demand  is  greater  than  1  [11]. 
Define  the  following.  .  . 


p  the  mean  of  the  sample  half-week  demands  divided  by  their  standard  package  sizes 
a2  — y  the  variance  of  the  sample  half-week  demands  divided  by  their  standard  package 
sizes 

p  — >  the  probability  of  success  approximated  as  p  =  -jfa 
r  — >  a  shape  factor  approximated  as  r  = 

k  — >■  the  number  of  failures  before  the  rth  success  in  a  sequence  of  independent  Bernoulli 
trials 


The  form  of  the  negative  binomial  [10]  that  was  used  follows. 


rrfc+V-r  —  n) 

P(Y  =  *)  =  ^  V(1  -  P)k  (5) 

Notice  that  r  is  allowed  to  take  on  real  values.  This  relaxed  version  allows  for  better  fits 
than  the  unrelaxed  distribution  and  is  usually  referred  to  as  a  generalized  negative  binomial 
distribution  [3]. 

Chi-square  goodness-of-fit  tests  unfortunately  do  not  unanimously  support  the  choice  of 
distribution.  Only  8  of  the  14  A  items  and  5  of  the  11  B  items  failed  to  reject  the  null 
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hypothesis  at  a  =  0.01  (critical  value  is  11.345  with  3  df)  (attachment  10).  One  logical 
explanation  is  the  erratic  nature  of  the  demand,  which  seems  to  prevent  it  from  being 
modeled  by  any  non-empirical  distribution.  Another  explanation  has  to  do  with  the  structure 
of  the  tests.  Using  6  intervals  and  having  to  estimate  2  distribution  parameters  leaves  only  3 
degrees  of  freedom  to  test  the  hypothesis.  Having  more  intervals  would  increase  the  degrees  of 
freedom  and  the  critical  chi-square  value,  possibly  resulting  in  some  rejection  reversals.  This, 
however,  leads  to  manipulating  the  tests  until  the  results  are  favorable,  a  time  consuming 
and  brow-raising  activity. 

The  other  initialization  requirement  for  this  heuristic  is  to  have  an  EOQ  value  to  work 
with  for  each  item.  Let  K  be  the  ordering  cost,  A  be  the  annual  demand  rate,  I  be  the 
holding  cost  rate,  and  c  be  the  standard  cost  for  some  item.  The  classic  EOQ  equation  then 
[1],  for  the  case  where  backorders  are  not  allowed,  is 


_  l2K\ 

E0Q=i~i -  (6) 

The  value  for  A  can  be  estimated  from  the  half-week  demands  in  the  sample.  The  c  values 
are  already  provided.  The  values  for  K  and  I  must  be  estimated  since  Quantegy  was  unable 
to  provide  such  information. 

Both  [1]  and  [8]  provide  equations  to  estimate  the  ratio  y  using  exchange  curves.  This 
approach  assumes  that  y  is  constant  across  all  items,  which  is  not  unreasonable.  It  also 
assumes  a  deterministic  system  where  all  items  are  being  replenished  with  their  EOQ  values. 
Obviously  this  is  not  happening  in  reality.  In  fact  most  of  the  EOQ  asumptions  (i.e.  no 
backorders,  infinite  production  rates,  constant  and  continuous  demand,  et  cetera)  are  inap¬ 
propriate  here.  But  the  EOQ  value  is  just  being  used  as  an  approximation  for  the  supposed 
best  value  for  S  —  s,  the  minimum  order  quantity  in  an  (s,  S )  inventory  control  system. 
And  making  the  EOQ  assumptions  incorporates  each  item’s  standard  cost,  something  that 
Quantegy’s  current  system  ignores.  This  heuristic  is  just  that:  a  heuristic.  It  does  not  claim 
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to  provide  optimal  s  and  S  values. 

The  exchange  curve  equation  that  relates  the  average  number  of  orders  per  year  over  all 
n  items  (call  it  R )  is  the  following  [1]  [8].  .  . 


To  estimate  R,  simulate  Quantegy’s  current  policy  over  the  sample  set  minus  the  first  13 
weeks  (used  to  initialize  the  demand  forecasts)  (attachment  11).  R  is  the  annualized  total 
number  of  orders  from  these  simlulations.  Substitute  this  value  into  the  equation  above  and 
solve  for  -j  (attachment  12).  The  resulting  value  of  y  is  19.535,  implying  that  the  ordering 
cost  varies  from  $1.95  to  $5.86  as  the  holding  cost  rate  varies  from  10%  to  30%.  EOQ  values 
can  now  be  calculated  for  every  item  (attachment  12). 

The  first  stage  can  now  begin.  Define  the  following  [14]  .  .  . 

•  pj  P [one-period  demand  =  j],  defined  by  the  distribution  above,  where  j  =  0, 1, 2... 

m(j)  the  expected  number  of  visits  to  state  y  -  j  until  the  next  order,  where  j  = 
0, 1,  ...y-  s  -  1 

M(j)  ->•  the  expected  number  of  periods  until  the  next  order  from  state  s  +  j,  where 
j  =  1  ,...S-s 

The  relevant  renewal  equations  for  the  first  stage  follow  below. 

j 

mti)  =  J2PimU~l)  (8) 

1=0 

M(j)  =  M(j  -  1)  +  m(j  -  1)  (9) 

where  m(0)  =  and  M( 0)  =  0. 

Recall  that  the  first  stage  selects  an  S  -  s  value  for  each  item  by  matching  its  EOQ  value 
to  the  corresponding  order  quantity  defined  by  the  renewal  equations.  This  latter  value  is 
defined  as  M(S  —  s)  *  p,  the  expected  cycle  length  (time)  multiplied  by  the  mean  demand 
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per  period  (quantity/time).  To  select  the  appropriate  S  -  s  value,  iterate  through  M(j), 
starting  at  M(  0),  until  M(j)  >  ^2.  The  resulting  j  value  is  S-s  for  that  item. 

Define  the  following  for  the  second  stage  [14]  ,  .  . 

A  -4,  the  S  -  s  value  found  in  the  first  stage 

L(y)  |  n  -¥  the  expected  number  of  stockouts  after  n  halfweeks  given  initial  inventory 
position  y 

H (s,  y)  -»■  the  expected  number  of  stockouts  until  the  next  order  given  initial  inventory 
position  y 

The  only  additional  renewal  equation  is  the  following. 

H  (s>  y)  =  £  mU)  [L(y  -j\2)~  L(y  -  j  I  1)]  (10) 

3=  0 

Representing  the  loss  function  in  the  equation  for  H(s,  y)  as  a  difference  between  the 
expected  loss  after  2  half-weeks  (lagtime  and  review  period)  and  the  expected  loss  after  1 
half-week  (lagtime  only)  avoids  the  common  mistake  of  double  counting  stockouts  when  the 
lagtime  is  not  taken  into  account.  Making  the  mistake  leads  to  underestimating  the  true 
fillrate.  This  was  pointed'out  by  Johnson,  et  al.  in  their  recent  paper  on  the  subject  [4], 

The  second  stage  of  this  heuristic  chooses  the  minimum  (s,  s  +  A)  policy  that  meets 
the  specified  fillrate.  Let  /?  be  the  specified  fillrate.  It  follows  that  the  expected  number  of 
stockouts  until  the  next  order  for  a  given  policy  that  achieves  the  fillrate  is  (1  -  /?)  *  M{ A)  *  y. 
To  select  the  appropriate  policy,  iterate  through  H(s,s  +  A),  starting  with  i/(0,  A)  until 
H(s,  s  +  A)  <  (1  —  (3)  *  M( A)  *  y.  The  resulting  policy  is  the  solution. 

The  (s,  S)  values  for  each  item  that  result  from  specifying  90%,  95%,  and  99%  fillrates 
are  given  in  attachment  13.  The  simulated  performance  of  these  static  policies  over  the  ex 
post  demands  follows  in  the  table  below  (attachment  14). 
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Performance  Measure 

P  =  90% 

0  =  95% 

P  =  99% 

Total  Orders 

997 

997 

997 

Avg  OH  Inv  ($) 

77,905.49 

106,028.00 

170,019.44 

#  Fillrates  >  90% 

18 

21 

25 

Avg  Fillrate 

92.7% 

96.2% 

99.1% 

Table  3:  Performance  Comparison  For  Two-Stage  Heuristic  at  Specified  Fillrates 

9  Alternative  3:  An  Almost-Exact  Algorithm 

The  section  details  an  algorithm  that  minimizes  an  ( s,S )  policy’s  average  ordering  and 

holding  costs  subject  to  a  minimum  expected  fillrate  constraint.  It  differs  from  the  problem 

Zheng  and  Federgruen  solved  in  [14]  in  that  here  the  backlogging  costs  associated  with  a 

policy  are  ignored,  essentially  being  replaced  with  the  fillrate  constraint. 

Define  the  following.  .  . 

h  -*•  the  annual  holding  cost  per  unit 
K  — v  the  fixed  cost  to  place  an  order 

G{y)  the  one-period  expected  holding  costs  for  a  given  policy  when  starting  from 
inventory  position  y 

c(s,  S )  — >  the  long-run  average  cost  for  a  given  policy 
P(s,  S )  ->  the  expected  fillrate  for  a  given  policy 
P  ->  the  minimum  acceptable  fillrate 

S  the  upperbound  on  S,  reflecting  the  storage  capacity  for  the  item  at  the  location 
ordering  it 

A  ->  the  upperbound  on  S  —  s,  reflecting  the  transportation  capacity  for  the  item  on  the 
vehicle  delivering  it 

The  relevant  equations  for  the  algorithm  follow  below. 


G{y)  =  h*[S  —  y  —  (2fj,  —  L(S  —  y)  |  2)] 

(11) 

c(s,  S)  =  M(S  -  s )~'[K  +  Sf;1  m(j)G(j)} 

3=0 

(12) 
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P(s,  S)  =  l  — 


H(s,S ) 


(13) 


M(S  -  s)v 

The  equation  for  G(y)  is  specific  for  Quantegy’s  case  where  the  lagtime  is  one  period  and 
the  review  interval  is  one  period.  It  captures  the  cost  of  holding  what  should  currently  be  on- 
hand  ( S-y )  minus  the  overshoot  caused  by  the  order  placed  a  week  prior  (2 y,-L(S-y)  |  2). 
The  equation  for  c(s,S)  can  apply  to  any  review  interval/lagtime  combination.  The  fillrate 
equation,  when  using  the  form  of  H(s,  y)  already  described,  is  only  applicable  for  Quantegy’s 
situation.  The  m(j)  and  M(j)  values  are  the  same  ones  used  in  the  prior  section. 

The  problem  then  is  to  minimize  c(s,  S )  subject  to 


.  i  0(8,  S)>P 

s<s 

S-s<  A 

The  following  lemmas  and  theorems  (proofs  in  attachment 

•  Lemma  1:  ( 3(s ,  S)  <  0(s  + 1,  S) 

•  Lemma  2:  ) 3(s ,  S )  <  0(s  +  1,  S  +  1) 

•  Lemma  3:  c(s,  S )  <  c(s  +  1,5) 

•  Lemma  4:  c(s,  S )  <  c(s  +  1,  S  +  1) 

•  Theorem  1:  If  (s*,  s*  +  1)  is  the  minimum  feasible  policy  such  that  /?(s*,  s*  +  1)  > 
then  all  policies  with  S  <  s*  are  infeasible. 

•  Theorem  2:  For  every  S  >  s*  +  1  there  is  one  s  for  which  (a)  all  policies  with  s  <  s 
and  S  <  S  are  infeasible  and  (b)  all  policies  with  s>  s  and  S'  >  5  are  dominated. 


The  algorithm  simply  searches  each  feasible  S  for  the  dominant  s  value  that  satisfies 
Theorem  2,  keeping  track  of  the  minimum  cost  policy  as  it  goes. 

\ 
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Exact  (s,  S)  Algorithm 


For  every  item.  .  . 

•  Search  for  the  (s*,  s*  +  1)  policy  defined  by  Theorem  1 

•  While  S  <  S 

—  While  [3(s,  S)  >  (3  and  s  >  S  —  A 
*  Decrement  s 

-  Increment  s  (become  feasible  again) 

-  If  c(s,S)  is  less  than  the  current  minimum  cost  policy,  then  (s,  S)  is  the  new 
current  minimum  cost  policy 

—  Increment  s  and  S 

Like  the  values  for  K  and  h,  the  upperbounds  S  and  A  are  unknown.  To  accomodate  this 
fact,  the  algorithm  is  modified  to  choose  the  first  policy  that  is  followed  by  20  consecutive 
non-improving  increments  of  S.  This  obviously  no  longer  results  in  an  optimal  solution, 
since  it  is  unknown  whether  policies  more  than  twenty  S  values  away  have  lower  cost  (hence 
the  “Almost- Exact” ) .  But  the  solutions  appear  to  be  optimal  anyway  because  using  a  value 
of  200  instead  of  20  did  not  result  in  any  difference  between  the  solutions.  [The  maximum 
number  of  states  examined  for  each  item  was  2000.] 

K  wasjrp.proxmmted  from  the  y  value  found  in  the  prior  section  by  assuming  a  holding 
costrateUilof  10%J  The  resulting  K  value  is  $1.95.  Using  the  same  I  value,  h  is  approx¬ 
imated  by  the  equation  h  =  Ic,  where  c  is  each  item’s  standard  cost.  It  turns  out  that  the 
algorithm  is  robust  with  respect  to  the  choice  of  /,  since  the  algorithm  produced  the  same 
policies  regardless  of  the  /  value  used  (non-zero  of  course).  The  specified  fillrates  again  are 
90%,  95%,  and  99%.  The  specific  policies  that  result  from  the  algorithm  are  included  in  at¬ 
tachment  16.  The  simulated  performance  of  the  policies  over  the  ex  post  set  are  summarized 
in  the  table  below  (attachment  17). 


Performance  Measure 

P  =  90% 

/?  =  95% 

P  =  99% 

Total  Orders 

1059 

1066 

1125 

Avg  OH  Inv  ($) 

77,141.52 

104,478.76 

168,584.11 

#  Fillrates  >  90% 

18 

21 

25 

Avg  Fillrate 

92.4% 

96.1% 

99.0% 

Table  4:  Performance  of  Almost-Exact  Alg.  at  Specified  Fillrates 


10  Recommendations  and  Conclusions 

Croston’s  method  was  the  clear  winner  of  the  forecasting  alternatives.  It  resulted  in  fewer 
orders  placed  and  superior  fillrates  than  Quantegy’s  current  method,  all  at  a  cost  of  a  mere 
$40  increase  in  on-hand  inventory.  The  naive  method  actually  produced  better  fillrates  than 
Croston’s  method,  but  had  unacceptably  high  levels  of  on-hand  inventory.  The  mean  forecast 
just  had  deplorable  fillrates.  The  resulting  conclusion  is  that  Quantegy  should  consider  using 
a  more  sophisticated  forecast  than  a  13-week  moving  average  with  their  current  system. 

The  renewal  analysis  approaches  at  yd  =  90%  both  produced  results  that  were  superior 
overall  to  Quantegy’s  current  system.  Both  the  two-stage  heuristic  and  almost-exact  poli¬ 
cies  had  fewer  orders,  lower  investments  in  on-hand  inventory,  and  higher  average  fillrates. 
Increasing  the  specified  fillrate  to  95%  and  99%  predictably  resulted  in  better  fillrate  perfor¬ 
mance,  but  pushed  inventory  costs  to  higher-than-current  levels.  As  a  result,  no  conclusions 
can  be  made  with  the  95%  and  99%  policies  because  of  the  subjective  nature  of  the  tradeoff 
between  fillrates  and  holding  costs.  The  90%  methods,  on  the  other  hand,  appear  to  be  clear 
recommendations  to  Quantegy’s  management,  with  the  almost-exact  policies  barely  edging 
out  their  two-stage  heuristic  competitors. 

There  is  one  drawback  to  these  recommendations  however.  Quantegy’s  current  simulated 
performance  has  all  14  of  the  A  items,  the  most  important,  meeting  at  least  a  90%  fillrate. 
Both  renewal  analysis  approaches,  on  the  other  hand,  have  only  8  of  the  14  A  items  meeting 
a  90%  fillrate.  This  may  be  unacceptable  to  Quantegy’s  service  goals  despite  the  overall 
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savings,  depending  on  just  how  important  the  A  items  are. 

The  first  intuitive  explanation  as  to  why  this  happens  deals  with  the  assumption  of  the 
underlying  distritbution.  Since  the  renewal  analysis  is  only  exact  if  the  specified  probability 
distribution  is  exact,  it  follows  that  the  items  that  failed  to  meet  90%  had  incorrect  distribu¬ 
tions.  The  chi-square  goodness-of-fit  (GOF)  tests,  however,  do  not  support  this  conclusion. 
Eight  of  the  14  A  items  “passed”  their  GOF  tests  with  respect  to  the  sample  demands.  Of 
these  same  8,  only  3  (A07,  A08,  A13)  still  “passed”  with  repsect  to  the  ex  post  demands 
(attachment  18).  Logically  these  3  items  should  have  met  the  90%  specified  fillrate.  But 
only  A08  did  so,  with  a  96.3%  fillrate. 

B  item  performance,  on  the  other  hand,  is  outstanding  for  the  renewal  analysis  ap¬ 
proaches.  At  j3  =  90%,  10  of  the.  11  B  items  meet  a  90%  fillrate  (B06).  This  is  despite  the 
fact  that  only  4  of  the  11  B  items  “pass”  the  GOF  test  for  the  ex  post  demands.  It  seems 
that  the  accuracies  of  the  simulated  fillrates  are  unpredictable  with  respect  to  the  GOF  tests. 

Perhaps  the  best  method  is  a  combination  of  both  Croston’s  method  and  the  almost-exact 
algorithm.  Specifically  the  almost-exact  algorithm  should  be  used  on  any  item  that  “passes” 
its  GOF  test  over  the  sample  set,  still  supporting  the  intuition  above.  The  remaining  items 
(A02-A05,  All,  A14,  B02,  B04-B06,  B08,  B10)  use  the  Croston  alternative.  The  table  below 
summarizes  the  simulated  performance  of  such  an  approach  (attachment  19).  Here  10  of  the 
14  A  items  now  meet  the  90%  requirement,  which  is  an  improvement  over  the  8  that  result 
from  renewal  analysis  alone.  The  investment  in  on-hand  inventory  is  the  second  lowest  so 
far.  The  number  of  orders  is  between  that  of  Croston’s  method  and  the  renewal  analysis 
approaches.  The  fillrate  performance  is  very  competitive.  Overall  this  approach  appears  to 
be  a  good  compromise  between  the  two  methods.  The  ultimate  decision  on  what  tradeoff  is 
best,  however,  and  what  approach  to  take  is  left  to  management. 

A  final  word  is  in  order  about  fairness.  The  renewal  analysis  policies  would  not  be 
static  in  reality.  The  analysts  at  Quantegy  would  most  likely  update  the  policies  monthly, 
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Performance  Measure 

Almost-Exact/Croston  Hybrid 

Total  Orders 

1100 

Avg  OH  Inv  ($) 

67,956.66 

#  Fillrates  >  90% 

17 

Avg  Fillrate 

92.0% 

Table  5:  Performance  of  Almost-Exact/Croston  Hybrid  Approach 


if  not  weekly,  as  new  information  came  in.  But  the  approach  in  this  paper  was  to  hold  the 
renewal  analysis  policies  constant  for  an  entire  year.  There  is  no  doubt  that  the  simulated 
performance  of  these  policies  would  have  improved  had  the  they  been  dynamic. 


11  Tasks  Before  Sumbitting  For  Publication 

The  extension  of  [14]  to  include  a  service  level  constraint  appears  to  be  non-trivial.  The 
main  reason  is  that  the  upperbound  on  S  that  the  authors  present  does  not  work  when 
using  a  service  level  constraint,  cutting  off  the  optimal  solution  too  soon  in  many  cases.  And 
the  upperbounds  presented  in  this  paper  (S, A)  are  questionable  both  in  their  existence  and 
availablilty.  For  example,  does  a  storage  limit  really  exist  for  any  one  individual  item  in  a 
warehouse?  If  it  did,  how  would  you  get  that  data? 

As  it  stands  now  there  are  no  formally  proven  upperbounds  /  stopping  criteria  for  the 
almost-exact  algorithm.  The  bound  with  the  most  promise  appears  to  be  the  minimum  (1,  S') 
policy  such  that  c(l,S)  >  c(s*,  s*  +  1),  where  s*  follows  the  notation  in  Theorem  1  above. 
Proving  that  c(l,  S  +  1)  >  c(l,  S)  under  these  assumptions  is  turning  out  to  be  harder  than 
expected.  Once  this  is  accomplished,  then  the  same  theorems  that  drive  the  algorithm  will 
stop  it,  resulting  in  an  algorithm  that  opt^ally  solves  the  (s,  S)  problem  subject  to  a  fillrate 
constraint. 


i  (ft 
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Raw  Data  -  A  Items 


Attachment  1 


API  A02  A03  A04  A05  A06  A07  A08  A09  A10  All  A12  A13  A14 

Mean  1723.90  290.90  61.83  53.08  127.21  125.14  412.62  355.52  146.29  247.28  130.27  39.80  191.10  480.81 

Standard  Dev  1121.93  260.55  65.28  56.26  120.24  85.37  219.34  231.24  155.91  177.08  138.18  38.81  164.17  814.15 

f  0.99  0.94  0.89  0.84  0.96  0.99  1.00  1.00  0.99  1.00  0.97  0.89  0.99  0.90 

Std  Package  10  10  5  5  10  10  10  10  10  1  2  7  7  10 
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B01  B02  B03  B04  BOS  B06  B07  B08 

Mean  81.56  1945  77.99  3^62  97^62  6Z23  29.48  163.60 

Standard  Dev  83.11  27.89  105.28  11.34  122.60  63.35  46.92  136.70 

f  0.93  0.69  0.84  0.30  0.88  0.86  0.67  0.97 
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Outlier 

Lower  Bound 
(Mean+3StdDev) 

Obs  >  Lower  Bound 

(Year  \  Week) 
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Value 

(Rounded  Mean) 
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5089.68 
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1720 

A02 
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257.66 

93\23,29  94\9 
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94\25,45  95\51 
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95\1  96\29,31 ,39 
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B01  B02  B03  B04  B05  B06  B07  B08  B09  BIO 

Mean  72.20  15.17  64.73  2.69  84.60  59.55  24.42  149.50  15.74  61.97 

Standard  Dev  61.55  17.31  65.42  5.43  80.70  52.73  31.79  96.56  32.89  97.98 

f  0.93  0.69  0.84  0.30  0.88  0.86  0.67  0.97  0.24  0.38 

Std  Package  10  5  10  1  10  10  10  2  40  20 
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A01 

A02 

A03 

A04 

A05 

A06 

A07 

A08 

A09 

A10 

All 

A12 

A13 

A14 

Sample  Mean 
Sample  S.D. 

762.11 

758.27 

118.40 

135.57 

32.56 

46.18 

28.64 

36.68 

47.89 

60.30 

52.86 

49.58 

200.42 

150.59 

174.21 

158.96 

71.41 

74.02 

126.35 

111.29 

70.93 

75.18 

18.43 

26.75 

78.29 

78.91 

182.60 

233.51 

Ex  Post  Mean 
Ex  Post  S.D. 

1009.41 

863.83 

174.61 

214.97 

18.82 

27.81 

13.92 

31.13 

80.39 

83.81 

78.43 

66.94 

198.24 

160.31 

157.98 

134.97 

44.12 

41.01 

89.28 

79.11 

29.61 

32.06 

20.38 

22.51 

95.19 

97.33 

198.73 

256.41 

Overall  Mean 
Overall  S.D. 

835.44 

797.82 

135.07 

164.87 

28.49 

42.02 

24.27 

35.72 

57.53 

69.60 

60.44 

56.42 

199.77 

153.30 

169.40 

152.23 

63.32 

67.09 

115.36 

104.08 

58.68 

68.06 

19.01 

25.55 

83.30 

84.99 

187.38 

240.26 

Half-Week  # 
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28.0  0  160  140  15  20  0 

28.5  260  120  45  85  10  100 

29.0  840  60  0  0  0  0 

29.5  170  440  0  25  90  40 

30.0  810  0  30  30  10  80 

30.5  920  0  20  60  50  0 

31.0  0  0  5  0  10  30 

31.5  2490  10  5  55  80  0 

32.0  1330  10  60  0  30  40 

32.5  470  0  5  0  0  0 

33.0  230  0  5  55  0  10 

33.5  670  10  60  45  170  10 

34.0  50  190  60  20  10  240 

34.5  360  280  75  0  0  10 

35.0  980  100  0  10'  10  20 

35.5  0  90  100  0  20  30 

36.0  380  10  15  75  40  20 

36.5  20  160  0  60  0  110 

37.0  3240  40  65  0  340  30 

37.5  1120  0  5  130  20  10 

38.0  0  280  5  0  0  20 

38.5  890  10  0  0  130  30 

39.0  150  250  60  0  0  30 

39.5  630  460  10  35  0  40 

40.0  470  0  75  0  90  60 

40.5  270  0  85  60  10  20 

41.0  490  100  65  0  30  100 

41.5  190  0  25  5  20  0 

42.0  30  60  0  45  50  10 

42.5  10  540  110  100  90  30 

43.0  870  40  50  0  60  40 

43.5  2980  0  50  55  0  10 

44.0  200  0  5  10  80  10 

44.5  640  200  0  0  0  40 

45.0  660  80  20  60  50  40 

45.5  1790  0  40  10  90  0 

46.0  30  0  0  110  20  60 

46.5  160  390  0  50  0  10 

47.0  540  270  140  20  30  60 

47.5  140  260  55  135  50  20 

48.0  420  0  30  35  0  20 

48.5  680  80  10  0  90  140 

49.0  680  40  20  5  30  50 

49.5  0  130  0  0  0  60 

50.0  600  70  10  60  40  50 

50.5  890  0  100  10  100  60 

51.0  1730  0  30  45  0  0 

51.5  680  0  20  60  0  20 

52.0  220  0  80  0  0  20 

52.5  320  140  65  20  50  0 

53.0  540  0  0  30  40  100 

53.5  720  90  20  0  110  80 

54.0  240  10  0  0  60  60 

54.5  0  10  0  35  10  90 

55.0  1400  290  21  20  30  10 

55.5  0  80  50  75  0  90 

56.0  620  20  0  0  80  0 

56.5  260  160  0  0  0  0 

57.0  2800  0  10  0  0  0 

57.5  370  250  35  0  20  20 

58.0  730  170  0  35  30  30 

58.5  400  20  15  55  80  10 

59.0  60  20  0  95  0  60 
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31.5  780  520  0  0  0  120 

32.0  2860  340  65  20  120  230 

32.5  1440  240  5  10  300  110 
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B01 

B02 

B03 

B04 

B05 

B06 

B07 

B08 

B09 

BIO 

B12 

Sample  Mean 

40.36 

9.67 

37.80 

1.38 

42.73 

23.57 

13.22 

82.25 

8.21 

34.71 

4.57 

Sample  S.D. 

52.39 

15.52 

57.94 

4.40 

59.44 

31.94 

24.84 

89.44 

24.73 

79.44 

15.66 

Ex  Post  Mean 

26.57 

2.65 

19.47 

1.26 

41.27 

44.51 

9.80 

56.96 

7.06 

22.16 

5.93 

Ex  Post  S.D. 

35.72 

5.48 

31.39 

2.97 

53.32 

59.15 

21.11 

60.51 

23.49 

56.84 

19.62 

Overall  Mean 

36.27 

7.59 

32.37 

1.34 

42.30 

29.78 

12.21 

74.75 

7.87 

30.99 

4.97 

Overall  S.D. 

48.41 

13.72 

52.15 

4.03 

57.62 

42.88 

23.81 

82.66 

24.34 

73.61 

16.92 

Half-Week  # 
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28.0  90  5  40  0  30 

28.5  30  0  0  0  40 

29.0  50  0  10  0  50 

29.5  60  50  30  0  120 

30.0  40  0  230  20  50 

30.5  40  10  40  0  10 

31.0  170  5  40  0  50 

31.5  30  0  0  0  230 

32.0  160  10  130  0  0 

32.5  0  5  0  0  30 

33.0  20  20  50  25  40 

33.5  10  0  20  0  0 

34.0  70  20  20  0  60 

34.5  10  30  40  0  20 

35.0  110  0  60  0  160 

35.5  0  0  10  0  120 

36.0  50  0  20  0  0 

36.5  40  20  80  0  0 

37.0  40  0  0  0  40 

37.5  40  0  20  0  0 

38.0  0  0  30  0  70 

38.5  0  0  0  0  240 

39.0  20  0  80  0  0 

39.5  0  45  0  0  20 

40.0  60  0  160  0  60 

40.5  0  0  20  0  100 

41.0  90  10  80  8  0 

41.5  200  0  40  0  20 

42.0  40  10  10  0  40 

42.5  40  0  10  0  40 

43.0  0  20  30  0  20 

43.5  70  0  140  0  20 

44.0  120  100  0  .0  70 

44.5  30  10  10  0  20 

45.0  50  ■  0  0  0  10 

45.5  10  0  0  0  0 

46.0  30  0  40  8  70 

46.5  0  5  0  0  0 

47.0  80  15  100  0  190 

47.5  10  5  180  0  0 

48.0  10  5  0  24  130 

48.5  160  0  0  0  60 
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50.5  30  15  160  0  40 
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54.5  0  0  220  0  84 

55.0  80  35  0  0  0 
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56.0  0  25  10  0  80 
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57.0  30  0  20  0  0 

57.5  10  35  350  0  0 

58.0  10  0  0  0  110 

58.5  30  20  80  0  40 

59.0  80  0  0  0  0 
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59.5  0  5  140  0  0  60  10  166  0  20  5 

60.0  240  0  0  0  180  0  0  34  0  140  15 

60.5  0  0  80  0  30  20  0  110  40  120  0 

61.0  160  20  30  0  20  20  0  90  0  0  0 

61.5  0  5  0  0  220  70  0  80  80  0  0 

62.0  90  5  100  0  0  0  0  114  0  0  20 

62.5  40  80  140  0  50  0  20  248  0  0  5 

63.0  30  0  10  0  20  20  0  14  0  0  5 

63.5  130  20  0  0  0  30  0  22  40  0  0 

64.0  0  0  60  0  0  70  0  108  0  200  0 

64.5  80  15  0  0  60  0  0  170  0  0  0 

65.0  0  0  250  0  0  50  10  108  80  160  0 

65.5  10  0  50  0  20  0  0  36  0  0  30 

66.0  0  0  30  0  0  10  10  66  0  0  0 

66.5  0  5  0  0  10  90  40  100  0  0  0 

67.0  50  0  160  0  20  0  50  64  0  240  20 

67.5  0  40  120  0  30  20  0  0  0  0  0 

68.0  130  20  0  0  0  20  0  144  0  0  0 

68.5  20  0  0  8  190  30  50  94  0  0  0 

69.0  0  0  100  0  100  0  0  152  0  0  0 

69.5  30  10  20  0  20  0  50  96  0  0  0 

70.0  40  0  30  0  40  110  20  0  0  0  5 

70.5  10  0  110  0  40  0  0  78  0  120  0 

71.0  10  0  10  0  30  20  0  26  0  0  0 

71.5  110  0  20  0  10  0  60  90  0  0  5 

72.0  60  35  90  0  0  40  0  34  40  320  0 

72.5  10  0  0  0  10  60  30  0  0  40  0 

73.0  70  35  10  0  40  0  0  280  0  120  120 

73.5  70  0  0  0  0  20  40  66  0  0  0 

74.0  60  5  70  0  20  10  0  208  80  0  0 

74.5  180  20  0  0  0  0  110  284  0  0  0 

75.0  10  0  60  0  0  0  0  32  40  0  0 

.75.5  60  10  0  0  0  20  0  170  0  0  0 

76.0  20  0  40  16  50  70  10  46  0  320  0 

76.5  0  5  0  0  100  90  0  160  0  0  0 

77.0  0  0  20  0  0  0  0  16  0  0  0 

77.5  10  15  10  0  60  20  30  214  0  0  0 

78.0  0  10  10  4  40  0  0  38  0  0  0 

78.5  0  0  10  8  0  30  10  0  40  0  0 

79.0  120  5  0  0  0  0  0  502  0  0  0 

79.5  20  60  40  0  50  0  10  164  0  120  0 

80.0  40  0  40  0  10  20  0  0  0  0  0 

80.5  40  25  30  0  110  0  0  2  0  0  0 

81.0  0  0  0  0  30  0  0  162  0  0  0 

81.5  60  20  10  0  80  0  20  0  40  0  20 

82.0  30  0  60  0  260  30  0  104  0  0  0 

82.5  90  5  60  0  110  30  40  50  0  0  0 

83.0  10  5  60  0  10  20  50  48  0  0  0 

83.5  50  25  0  0  10  20  20  22  0  0  0 

84.0  0  0  0  0  20  90  0  100  0  0  0 

84.5  80  0  20  0  60  0  0  0  108  0  0 

85.0  0  5  0  0  60  20  0  118  0  0  0 

85.5  0  0  0  0  0  40  0  126  0  0  0 

86.0  20  40  110  0  90  0  0  0  0  0  0 

86.5  30  10  0  0  230  20  30  66  0  0  0 

87.0  170  15  230  0  0  50  0  26  0  0  15 

87.5  40  50  30  0  40  0  0  22  0  0  0 

88.0  0  0  0  10  100  40  20  182  0  0  0 

88.5  0  10  0  0  10  0  0  74  0  0  20 

89.0  10  15  0  24  10  0  0  56  0  0  0 

89.5  220  25  0  0  50  30  100  10  0  0  90 

90.0  70  0  280  24  60  0  0  172  120  0  0 

90.5  0  20  10  4  30  70  0  108  0  0  0 
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0 

90 

0 

110 

160 

20 

46 

0 

0 

95 

48.0 

20 

0 

40 

0 

0 

0 

0 

148 

0 

0 

0 

48.5 

0 

0 

30 

0 

20 

170 

0 

172 

0 

0 

0 

.49.0 

0 

0 

10 

0 

20 

30 

0 

0 

40 

0 

0 

49.5 

40 

0 

70 

0 

160 

0 

70 

0 

0 

0 

0 

50.0 

0 

0 

10 

15 

120 

40 

0 

142 

0 

0 

0 

50.5 

0 

15 

0 

0 

0 

10 

0 

0 

0 

0 

0 

51.0 

0 

0 

0 

0 

20 

0 

0 

82 

0 

0 

0 

Simulation  Results  -  Existing  Policy 


Attachment  6 


c 

Avg  OH  Inv 

$  OH  Inv 

#  Orders 

Fillrate 

Meet  90%? 

A01 

5.475 

4142.549 

22680.46 

61 

1.000 

1 

A02 

7.075 

665.882 

4711.12 

53 

0.964 

1 

A03 

31.205 

71.373 

2227.19 

48 

0.904 

1 

A04 

33.906 

83.284 

2823.83 

24 

0.954 

1 

A05 

6.049 

270.980 

1639.16 

63 

0.996 

1 

A06 

6.703 

283.235 

1898.52 

66 

0.998 

1 

A07 

9.161 

793.333 

7267.72 

70 

1.000 

1 

A08 

17.468 

682.275 

11917.98 

66 

1.000 

1 

A09 

10.987 

161.275 

1771.93 

65 

0.998 

1 

A10 

17.164 

345.745 

5934.37 

65 

0.991 

1 

All 

57.509 

101.529 

5838.83 

61 

0.942 

1 

A12 

35.970 

80.637 

2900.51 

52 

0.997 

1 

A13 

20.977 

330.922 

6941.75 

65 

0.952 

1 

A14 

5.292 

775.098 

4101.82 

50 

0.974 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  1 

82655.19 

809 

0.976 

1.000 

B01 

9.744 

94.216 

918.04 

46 

0.830 

0 

B02 

22.467 

10.490 

235.68 

22 

0.778 

0 

B03 

7.689 

62.882 

483.50 

40 

0.785 

0 

B04 

63.440 

4.706 

298.55 

25 

0.636 

0 

B05 

10.928 

121.177 

1324.22 

51 

0.843 

0 

B06 

11.360 

116.863 

1327.56 

52 

0.879 

0 

B07 

6.916 

35.980 

248.84 

33 

0.750 

0 

B08 

70.470 

140.569 

9905.90 

68 

0.935 

1 

B09 

3.106 

31.373 

97.44 

11 

0.556 

0 

BIO 

3.576 

86.471 

309.22 

23 

0.602 

0 

B12 

28.399 

21.275 

604.19 

20 

0.430 

0 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  1 

15753.14 

391 

0.729 

0.091 

Total 

Total 

Average 

Total  %  Meeting  90% 

98408.33 

1200 

0.868 

0.600 

Initialized  Croston's  Method  Parameters 


Attachment  7 


Alpha 

Z 

d 

P 

A01 

0.053 

1942.13 

1 

1.0003 

A02 

0.035 

276.06 

1 

1.0205 

A03 

0.090 

45.64 

2 

1.1597 

A04 

0.060 

62.64 

2 

1.0740 

A05 

0.093 

110.92 

1 

1.0304 

A06 

0.078 

138.40 

1 

1.0005 

A07 

0.047 

442.12 

1 

1.0000 

A08 

0.172 

437.08 

1 

1.0000 

A09 

0.074 

78.26 

1 

1.0006 

A10 

0.049 

201.34 

1 

1.0000 

All 

0.141 

41.36 

1 

1.1690 

A12 

0.032 

39.25 

1 

1.1136 

A13 

0.014 

153.71 

1 

1.0251 

A14 

0.063 

403.31 

1 

1.0651 

B01 

0.013 

84.53 

1 

1.0360 

B02 

0.053 

22.14 

2 

1.5047 

B03 

0.028 

79.49 

1 

1.1509 

B04 

0.076 

9.11 

4 

3.0303 

B05 

0.008 

93.19 

1 

1.1648 

B06 

0.047 

53.57 

1 

1.1521 

B07 

0.030 

36.88 

1 

1.4011 

B08 

0.164 

79.86 

1 

1.1776 

B09 

0.143 

56.91 

1 

3.4576 

BIO 

0.053 

165.17 

3 

2.5910 

B12 

0.026 

30.27 

3 

2.6915 

Forecast  Accuracies 


Attachment  8 


MvAvg  MSE 


1121244.51 

70674.52 

1288.07 

1703.79 

8090.67 

6220.22 

49696.75 

29552.67 
3319.97 

11705.68 
1829.78 
1040.44 

15015.78 

79470.31 


2425.25 

61.18 

2035.08 

17.39 

6490.78 

4332.85 

864.66 

5016.01 

1071.49 

6184.34 

769.61 


Overall  Average  %  Difference 


Naive  MSE 

%  Diff 

Mean  MSE 

%  Diff 

2189474.51 

0.95 

2622322.19 

1.34 

172766.67 

1.44 

120198.47 

0.70 

2278.92 

0.77 

1207.58 

-0.06 

3498.04 

1.05 

1738.13 

0.02 

15537.25 

0.92 

22034.13 

1.72 

10888.24 

0.75 

16621.66 

1.67 

98488.24 

0.98 

83616.57 

0.68 

50387.84 

0.71 

47289.51 

0.60 

6507.84 

0.96 

3211.24 

-0.03 

20807.71 

0.78 

13351.72 

0.14 

2906.90 

0.59 

1672.23 

-0.09 

2418.29 

1.32 

1463.10 

0.41 

23305.75 

0.55 

26278.86 

0.75 

189066.67 

1.38 

117962.83 

0.48 

Avg  %  Diff 

0.94 

Avg  %  Diff 

0.60 

3805.88 

0.57 

2290.58 

-0.06 

131.86 

1.16 

76.41 

0.25 

3735.92 

0.84 

1868.33 

-0.08 

27.33 

0.57 

16.52 

-0.05 

14476.47 

1.23 

7408.36 

0.14 

7307.84 

0.69 

8496.92 

0.96 

1915.69 

1.22 

840.61 

-0.03 

9250.82 

0.84 

6267.06 

0.25 

2227.45 

1.08 

1027.69 

-0.04 

11278.43 

0.82 

5697.22 

-0.08 

1612.75 

1.10 

745.36 

-0.03 

Avg  %  Diff 

0.92 

Avg  %  Diff 

0.11 

nee 

0.93 

0.38 

Croston  MSE 


1084032.66 

70187.93 

1261.72 

1799.73 
8529.13 
6065.49 

47345.25 
26922.41 
3117.15 
10775.47 
1659.12 
1000.54 
14720.43 
75355.84 


Avg  %  Diff 


2612.46 


2218.13 


5863.66 

4462.99 

834.26 

4891.71 

1012.43 

5856.01 


Avg  %  Diff 


%  Diff 


-0. 

03 

-0. 

01 

-0. 

02 

0. 

06 

0. 

05 

-0. 

02 

-0. 

05 

-0. 

09 

-0. 

06 

-0. 

08 

-0. 

09 

Simulation  Results  -  Naive  Forecast 


Attachment  9 


Total  Total  Average  Total  %  Meeting  90% 
173144.71  590  0.958  0.880 


Simulation  Results  -  Mean  Forecast 


Attachment  9 


c 

Avg  OH  Inv 

$  OH  Inv 

#  Orders 

Fillrate 

Meet  90%? 

A01 

5.475 

620.980 

3399.87 

82 

0.693 

0 

7.075 

120.882 

855.24 

69 

0.562 

0 

31.205 

61.814 

1928.91 

42 

0.924 

1 

A04 

33.906 

59.853 

2029.38 

27 

0.799 

0 

A05 

6.049 

31.373 

189.78 

78 

0.521 

0 

A06 

6.703 

36.177 

242.49 

79 

0.660 

0 

A07 

9.161 

230.490 

2111.52 

88 

0.868 

0 

A08 

17.468 

214.333 

3743.97 

78 

0.932 

1 

A09 

10.987 

119.216 

1309.83 

68 

0.987 

1 

A10 

17.164 

201.735 

3462.58 

72 

0.921 

1 

All 

57.509 

148.961 

8566.60 

53 

1.000 

1 

A12 

35.970 

20.931 

752.89 

57 

0.724 

0 

A13 

20.977 

79.265 

1662.74 

77 

0.737 

0 

A14 

5.292 

211.471 

1119.10 

52 

0.794 

0 

Subtotal  Subtotal  Subaverage  %  Meeting  90% 
31374.88  922  0.794  0.357 

B01  9.744  53.922  525.42  51  0.731  0 


B02  22.467  18.529  416.29  17  0.889  0 

B03  7.689  56.804  436.77  53  0.780  0 

B04  63.440  1.824  115.71  23  0.403  0 

B05  10.928  44.804  489.62  59  0.672  0 

B06  11.360  15.882  180.42  53  0.359  0 

B07  6.916  17.353  120.01  37  0.600  0 

B08  70.470  98.333  6929.53  69  0.903  1 

B09  3.106  15.686  48.72  6  0.389  0 

BIO  3.576  55.294  197.73  20  0.496  0 

B12 _ 28.399 _ 7.647  217.17 _ 13  0.198 _ 0 _ 

Subtotal  Subtotal  Subaverage  %  Meeting  90% 

_ I _ 9677.38  401  0.584 _ 0.091 

Total  Total  Average  Total  %  Meeting  90% 
41052.27  1323  0.702  0.240 


Simulation  Results  -  Croston  Forecast 


Attachment  9 


Subtotal 


Subtotal 


Subtotal 

16946.51 


Total 

98448.10 


Fillrate  Meet  90%? 


1.000  1 

0.957  1 

0.930  1 

0.975  1 

0.994  1 

1.000  1 

1.000  1 

1.000  1 

1.000  1 

1.000  1 

0.974  1 

0.983  1 

0.946  1 

0.984  1 


Subaverage  %  Meeting  90% 

0.982  1.000 


1 

l  1 

1 

i  0 


Subtotal 

358 


Total 

1151 


0.941 

1 

0.611 

0 

0.761 

0 

0.463 

0 

Subaverage  %  Meeting  90% 

0.809 

0.364 

Average 

Total  %  Meeting  90% 

0.906 

0.720 

Sample  Set  GOF  Test  Results  -  A  items 


Attachment  10 


A01 


Interval 

Prob 

0 

E 

Stat 

Othr  13 

0.1620 

38 

39.2120 

0.0375 

14  thr  30 

0.1656 

36 

40.0833 

0.4160 

31  thr  52 

0.1680 

46 

40.6573 

0.7021 

53  thr  82 

0.1643 

38 

39.7634 

0.0782 

83  thr  133 

0.1666 

41 

40.3057 

0.0120 

134  thr  inf 

0.1735 

43 

41.9783 

0.0249 

1.2705 

A02 

Interval 

Prob 

0 

E 

Stat 

0  thr  1 

0.1884 

77 

45.6013 

21.6195 

2  thr  3 

0.1300 

18 

31.4704 

5.7658 

4  thr  7 

0.1878 

25 

45.4587 

9.2075 

8  thr  12 

0.1567 

25 

37.9126 

4.3979 

13  thr  21 

0.1630 

52 

39.4515 

3.9914 

22  thr  inf 

0.1740 

45 

42.1056 

0.1990 

45.1810 

A03 

Interval 

Prob 

0 

E 

Stat 

0 

0.2504 

83 

60.5853 

8.2928 

1 

0.1245 

13 

30.1271 

9.7367 

2  thr  3 

0.1576 

34 

38.1316 

0.4477 

4  thr  6 

0.1442 

29 

34.8945 

0.9957 

7  thr  12 

0.1542 

42 

37.3152 

0.5882 

13  thr  inf 

0.1692 

41 

40.9462 

0.0001 

20.0611 

A04 

Interval 

Prob 

0 

E 

Stat 

0 

0.2170 

98 

52.5087 

39.4118 

1 

0.1322 

16 

31.9990 

7.9993 

2  thr  3 

0.1788 

18 

43.2578 

14.7478 

4  thr  6 

0.1662 

22 

40.2127 

8.2487 

7  thr  11 

0.1496 

41 

36.2011 

0.6361 

12  thr  inf 

0.1563 

47 

37.8207 

2.2279 

73.2716 

AOS 

Interval 

Prob 

0 

E 

Stat 

0 

0.2293 

71 

55.4965 

4.3311 

1 

0.1446 

32 

35.0045 

0.2579 

2 

0.1084 

12 

26.2368 

7.7253 

3  thr  4 

0.1548 

32 

37.4515 

0.7935 

5  thr  9 

0.2038 

51 

49.3240 

0.0569 

10  thr  inf 

0.1590 

44 

38.4867 

0.7898 

13.9545 

A06 

Interval 

Prob 

0 

E 

Stat 

0  thr  1 

0.2307 

60 

55.8387 

0.3101 

2 

0.1179 

30 

28.5436 

0.0743 

3  thr  4 

0.1993 

38 

48.2329 

2.1710 

5  thr  6 

0.1465 

36 

35.4481 

0.0086 

7  thr  9 

0.1406 

37 

34.0320 

0.2588 

10  thr  inf 

0.1649 

41 

39.9046 

0.0301 

2.8529 

A07 

Interval 

Prob 

0 

E 

Stat 

Othr  6 

0.1645 

48 

39.8133 

1.6834 

7  thr  11 

0.1729 

41 

41.8367 

0.0167 

12  thr  16 

0.1607 

30 

38.8837 

2.0297 

17  thr  23 

0.1780 

32 

43.0854 

2.8522 

24  thr  32 

0.1493 

44 

36.1229 

1.7177 

33  thr  inf 

0.1746 

47 

42.2580 

0.5321 

8.8318 

A08 


Interval 

Prob 

0 

E 

Stat 

Othr  3 

0.1520 

41 

36.7842 

0.4832 

4  thr  8 

0.1980 

48 

47.9185 

0.0001 

9  thr  13 

0.1640 

36 

39.6835 

0.3419 

14  thr  19 

0.1490 

34 

36.0677 

0.1185 

20  thr  30 

0.1695 

37 

41.0140 

0.3928 

31  thr  inf 

0.1675 

46 

40.5321 

0.7376 

2.0742 

A09 

Interval 

Prob 

0 

E 

Stat 

0 

0.1130 

31 

27.3374 

0.4907 

1  thr  2 

0.1998 

43 

48.3454 

0.5910 

3  thr  5 

0.2245 

58 

54.3296 

0.2480 

6  thr  8 

0.1532 

34 

37.0758 

0.2552 

9  thr  13 

0.1523 

36. 

36.8581 

0.0200 

14  thr  inf 

0.1572 

40 

38.0537 

0.0995 

1.7044 

A10 

Interval 

Prob 

0 

E 

Stat 

0  thr  30 

0.1623 

45 

39.2657 

0.8374 

31  thr  59 

0.1642 

38 

39.7402 

0.0762 

60  thr  95 

0.1726 

37 

41.7654 

0.5437 

96  thr  138 

0.1561 

34 

37.7779 

0.3778 

139  thr  21 5 

0.1733 

42 

41.9442 

0.0001 

216  thr  inf 

0.1715 

46 

41.5066 

0.4864 

2.3217 

All 

Interval 

Prob 

0 

E 

Stat 

Othr  5 

0.1720 

67 

41.6143 

15.4858 

6  thr  12 

0.1503 

18 

36.3752 

9.2824 

13  thr  23 

0.1775 

36 

42.9664 

1.1295 

24  thr  38 

0.1665 

29 

40.2926 

3.1649 

39  thr  63 

0.1618 

48 

39.1552 

1.9979 

64  thr  inf 

0.1719 

44 

41.5962 

0.1389 

31.1995 

A12 

Interval 

Prob 

0 

E 

Stat 

0 

0.3710 

98 

89.7834 

0.7520 

1 

0.1760 

38 

42.5878 

0.4942 

2 

0.1139 

26 

27.5572 

0.0880 

3  thr  4 

0.1391 

29 

33.6598 

0.6451 

5  thr  7 

0.1038 

24 

25.1190 

0.0499 

8  thr  inf 

0.0963 

27 

23.2929 

0.5900 

2.6191 

A13 

Interval 

Prob 

0 

E 

Stat 

Othr  1 

0.1440 

45 

34.8511 

2.9554 

2  thr  4 

0.1901 

36 

46.0102 

2.1779 

5  thr  8 

0.1946 

40 

47.0906 

1.0676 

9  thr  12 

0.1396 

34 

33.7933 

0.0013 

13  thr  20 

0.1689 

48 

40.8623 

1.2468 

21  thr  inf 

0.1628 

39 

39.3925 

0.0039 

7.4529 

A14 

Interval 

Prob 

0 

E 

Stat 

0 

0.1166 

83 

28.2161 

106.3672 

1  thr  4 

0.2172 

26 

52.5532 

13.4164 

5  thr  10 

0.1829 

29 

44.2619 

5.2625 

11  thr  17 

0.1358 

11 

32.8685 

14.5498 

18  thr  34 

0.1795 

37 

43.4342 

0.9531 

35  thr  inf 

0.1680 

56 

40.6660 

5.7820 

146.3310 


Sample  Set  GOF  Test  Results  -  B  Items 
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B01 


Interval 

Prob 

O 

E 

Stat 

0 

0.2634 

77 

63.7385 

2.7592 

1 

0.1564 

25 

37.8391 

4.3564 

2 

0.1131 

24 

27.3687 

0.4146 

3  thr  4 

0.1550 

42 

37.5112 

0.5372 

5  thr  8 

0.1649 

41 

39.8965 

0.0305 

9  thr  inf 

0.1473 

33 

35.6460 

0.1964 

8.2943 

B02 

Interval 

Prob 

0 

E 

Stat 

0 

0.4584 

129 

110.9239 

2.9457 

1 

0.1780 

27 

43.0763 

5.9998 

2 

0.1057 

18 

25.5773 

2.2448 

3 

0.0700 

13 

16.9383 

0.9157 

4  thr  5 

0.0837 

28 

20.2560 

2.9606 

6  thr  inf 

0.1042 

27 

25.2282 

0.1244 

15.1909 

B03 

Interval 

Prob 

O 

E 

Stat 

0 

102 

84.8960 

1 

0.1493 

29 

36.1349 

1.4088 

2 

0.0980 

15 

23.7233 

3.2076 

3  thr  4 

0.1274 

29 

30.8331 

0.1090 

5  thr  8 

0.1336 

31 

32.3431 

0.0558 

9  thr  inf 

0.1408 

36 

34.0697 

0.1094 

8.3365 

B04 

Interval 

Prob 

0 

E 

Stat 

0 

0.7573 

209 

183.2654 

3.6137 

1 

0.0739 

1 

17.8893 

15.9452 

2 

0.0379 

0 

9.1833 

9.1833 

3  thr  4 

0.0426 

11 

10.3043 

0.0470 

5  thr  9 

0.0463 

10 

11.2030 

0.1292 

10  thr  inf 

0.0420 

11 

10.1548 

0.0704 

28.9887 

B05 

Interval 

Prob 

0 

E 

Stat 

0 

91 

69.8963 

6.3718 

1 

0.1493 

19 

36.1230 

8.1166 

2 

0.1042 

23 

25.2114 

0.1940 

3  thr  4 

0.1413 

29 

34.1924 

0.7885 

5  thr  8 

0.1538 

39 

37.2154 

0.0856 

9  thr inf 

0.1627 

41 

39.3615 

0.0682 

15.6247 

B06 

Interval 

Prob 

0 

E 

Stat 

0 

0.3544 

108 

85.7678 

5.7629 

1 

0.1929 

20 

46.6884 

15.2558 

2 

0.1267 

31 

30.6595 

0.0038 

3 

0.0879 

24 

21.2815 

0.3473 

4  thr  5 

0.1081 

18 

26.1550 

2.5427 

6  thr  inf 

0.1299 

41 

31.4478 

2.9015 

26.8140 


B07 


Interval 

Prob 

0 

E 

Stat 

0 

0.5737 

147 

138.8417 

0.4794 

1 

0.1626 

34 

39.3492 

0.7272 

2 

0.0869 

15 

21.0338 

1.7309 

3 

0.0537 

12 

13.0041 

0.0775 

4  thr  5 

0.0598 

20 

14.4661 

2.1170 

6  thr  inf 

0.0632 

14 

15.3050 

0.1113 

5.2432 

B08 

Interval 

Prob 

0 

E 

Stat 

0  thr  5 

0.1633 

59 

39.5181 

6  thr  14 

0.1701 

31 

41.1584 

2.5072 

15  thr  26 

0.1654 

29 

40.0387 

3.0434 

27  thr  44 

0.1693 

28 

40.9791 

4.1108 

45  thr  74 

0.1617 

47 

39.1381 

1.5793 

75  thr  inf 

0.1701 

48 

41.1676 

1.1339 

21.9790 

B09 

Interval 

Prob 

0 

E 

Stat 

0 

211 

208.6716 

1 

0.0952 

18 

23.0347 

1.1004 

2 

0.0273 

9 

6.5981 

0.8743 

3 

0.0094 

2 

2.2772 

0.0337 

4 

0.0035 

2 

0.8528 

1.5434 

5  thr  inf 

0.0023 

0 

0.5656 

0.5656 

4.1435 

BIO 

Interval 

Prob 

0 

E 

Stat 

0 

0.6229 

187 

150.7400 

8.7222 

1 

0.1189 

2 

28.7712 

24.9102 

2 

0.0643 

8 

15.5491 

3.6651 

3  thr  4 

0.0724 

8 

17.5213 

5.1740 

5  thr  7 

0.0540 

10 

13.0701 

0.7212 

8  thr  inf 

0.0676 

27 

16.3483 

mm 

6.9402 

50.1328 

B12 

Interval 

Prob 

0 

E  ’ 

Stat 

0 

0.8005 

202 

193.7263 

1 

0.0680 

9 

16.4625 

3.3828 

2 

0.0337 

5 

8.1648 

1.2267 

3  thr  4 

0.0363 

13 

8.7930 

2.0128 

5  thr  8 

0.0319 

6 

7.7167 

0.3819 

9  thr  inf 

0.0295 

7 

7.1367 

0.0026 

7.3601 


Simulation  Results  -  Existing  Policy  Over  Sample  Set  (minus  first  1 3  weeks)  Attachment  1 1 


c 

Avg  OH  Inv 

$  OH  Inv 

#  Orders 

Fillrate 

Meet  90%  ? 

A01 

5.475 

2886.296 

15802.47 

137 

0.990 

1 

7.075 

446.398 

3158.27 

114 

0.968 

1 

EHi 

31.205 

139.579 

4355.55 

96 

0.989 

1 

A04 

33.906 

115.093 

3902.33 

104 

0.974 

1 

A05 

6.049 

195.417 

1182.08 

113 

0.965 

1 

A06 

6.703 

201.870 

1353.14 

133 

0.983 

1 

A07 

9.161 

774.556 

7095.70 

143 

1.000 

1 

A08 

17.468 

672.593 

11748.85 

143 

0.997 

1 

A  09 

10.987 

299.426 

3289.79 

120 

0.985 

1 

A10 

17.164 

519.954 

8924.49 

143 

0.986 

1 

All 

57.509 

326.296 

18764.97 

112 

0.997 

1 

A12 

35.970 

75.509 

2716.07 

93 

0.930 

1 

A13 

20.977 

315.208 

6612.12 

130 

0.980 

1 

A14 

5.292 

724.398 

3833.51 

106 

0.972 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  j 

92739.34 

1687 

0.980 

1.000 

B01 

9.744 

127.630 

1243.62 

106 

0.922 

1 

22.467 

32.917 

739.54 

80 

0.872 

0 

B03 

7.689 

129.444 

995.30 

92 

0.853 

0 

B04 

63.440 

'  7.005 

444.37 

37 

0.511 

0 

B05 

10.928 

139.769 

1527.39 

104 

0.839 

0 

B06 

11.360 

70.194 

797.41 

1 0S- 

0.928 

1 

B07 

6.916 

43.102 

298.09 

72 

0.799 

0 

B08 

70.470 

272.444 

19199.16 

122 

0.981 

1 

B09 

3.106 

34.815 

108.13 

26 

0.591 

0 

BIO 

3.576 

142.875 

510.92 

57 

0.759 

0 

B12 

28.399 

27.917 

792.81 

31 

0.565 

0 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  I 

26656.74 

832 

0.784 

0.273 

Total 

Total 

Average 

Total  %  Meeting  90% 

• 

119396.09 

2519 

0.893 

0.680 

Exchange  Curve  /  EOQ  Calculations 


Attachment  12 


c 

Mean 

Hfwk  Demand 

lambda 

SQRT  of 
c*lambda 

Resulting 
EOQ  Value 

A01 

5.475 

775.694 

79120.829 

658.169 

751.407 

A02 

7.075 

119.602 

12199.394 

293.787 

259.554 

A03 

31.205 

29.144 

2972.639 

304.567 

61.007 

A04 

33.906 

29.421 

3000.973 

318.984 

58.805 

A05 

6.049 

50.231 

5123.611 

176.048 

181.915 

A06 

6.703 

54.551 

5564.195 

193.124 

180.090 

A07 

9.161 

205.653 

20976.586 

438.368 

299.101 

A08 

17.468 

187.732 

19148.613 

578.349 

206.952 

A09 

10.987 

68.991 

7037.055 

278.058 

158.190 

A10 

17.164 

125.435 

12794.390 

468.618 

170.656 

All 

57.509 

66.333 

6766.000 

623.784 

67.799 

A12 

35.970 

17.986 

1834.583 

256.885 

44.640 

A13 

20.977 

78.315 

7988.111 

409.349 

121.975 

A14 

5.292 

191.250 

19507.500 

321.300 

379.501 

B01 

9.744 

39.713 

4050.722 

198.671 

127.444 

B02 

22.467 

9.051 

923.194 

144.019 

40.068 

B03 

7.689 

37.546 

3829.723 

171.601 

139.499 

B04 

63.440 

1.431 

145.917 

96.213 

9.480 

B05 

10.928 

43.764 

4463.917 

220.866 

126.331 

B06 

11.360 

22.514 

2296.417 

161.516 

88.871 

B07 

6.916 

12.917 

1317.500 

95.456 

86.272 

B08 

70.470 

79.648 

8124.111 

756.641 

67.113 

B09 

3.106 

8.463 

863.222 

51.780 

104.204 

BIO 

3.576 

34.995  • 

3569.528 

112.981 

197.483 

B12 

28.399 

3.889 

396.667 

106.136 

23.361 

Sum  of  SQRTsj5257.529  | 

Total  #  orders  from  Sample  Set  sim  2519 

Avg  #  Replenishments  /  Year  1 189.5278 
K/l  EstimatefT9~535034  | 


Two-Stage  Heuristic  Policies 


Attachment  13 


Simulation  Results  -  Two-Stage  Heuristic  at  90% 


Attachment  14 


c 

Avg  OH  Inv 

$  OH  Inv 

5.475 

1158.137 

6340.80 

7.075 

250.098 

1769.44 

31.205 

115.147 

3593.17 

33.906 

114.902 

3895.87 

6.049 

111.569 

674.88 

6.703 

82.059 

550.04 

9.161 

279.412 

2559.69 

17.468 

349.824 

6110.72 

10.987 

193.628 

2127.39 

17.164 

261.490 

4488.22 

57.509 

214.824 

12354.28 

35.970 

54.284 

1952.61 

20.977 

135.882 

2850.41 

5.292 

542.255 

2869.61 

Subtotal 

52137.12 

9.744 

146.275 

1425.30 

22.467 

50.245 

1128.86 

7.689 

183.471 

1410.71 

63.440 

'21.324 

1352.76 

10.928 

150.000 

1639.20 

11.360 

55.392 

629.25 

6.916 

79.118 

547.18 

70.470 

197.843 

13942.00 

3.106 

117.255 

364.19 

3.576 

308.039 

1101.55 

28.399 

78.431 

2227.37 

Subtotal 

25768.38 

#  Orders 


30 

19 

38 

41 

63 
83 
32 

64 
79 
44 

65 

42 


Subtotal 

747 


21 


Fillrate 


0.861 

0.759 

0.990 

0.958 

0.793 

0.809 

0.897 

0.965 

1.000 

0.946 

1.000 

0.953 

0.851 

0.962 


Meet  90%? 


0 

0 

1 

1 

0 

0 

0 

1 

1 

1 

1 

1 

0 

1 


Total 

77905.49 


14 

11 


Subtotal 

250 


Total 

997 


Subaverage  %  Meeting  90% 


Gian 


1.000  1 

1.000  1 

1.000  1 

0.922  1 

0.711  0 

0.940  1 

0.983  1 

1.000  1 

0.973  1 

0.942  1 


Subaverage  %  Meeting  90% 

0.949  0.909 


Average  Total  %  Meeting  90% 
0.927  0.720 


Simulation  Results  -  Two-Stage  Heuristic  at  95% 


Attachment  14 


5.475 

7.075 

31.205 

33.906 

6.049 

6.703 

9.161 

17.468 

10.987 

17.164 

57.509 

35.970 

20.977 

5.292 


9.744 

22.467 

7.689 

63.440 

10.928 

11.360 

6.916 

70.470 

3.106 

3.576 

28.399 


Avg  OH  Inv 


1724.900 

345.784 

154.951 

148.775 

153.333 

115.392 

387.941 

476.490 

263.627 

345.931 

278.824 

81.186 

197.578 

756.275 


195.294 

65.245 

253.471 

29.324 

208.039 

69.706 

108.529 

264.980 

157.255 

427.451 

107.745 


$  OH  Inv 


9443.83 

2446.42 

4835.25 

5044.37 

927.51 

773.47 

3553.93 

8323.33 

2896.47 

5937.56 

16034.89 

2920.27 

4144.59 

4002.21 


Subtotal 

71284.09 


1902.94 
1465.86 

1948.94 
1860.28 
2273.45 

791.86 

750.59 

18673.14 

488.43 

1528.56 

3059.85 


Subtotal 

34743.91 


Total 

106028.00 


#  Orders 


30 

19 

38 

41 

63 
83 
32 

64 
79 
44 

65 

42 


Subtotal 

747 


14 

11 


Subtotal 

250 


Total 

997 


Fillrate 

Meet  90%? 

0.945 

1 

0.860 

0 

1.000 

1 

1.000 

1 

0.871 

0 

0.876 

0 

0.941 

1 

0.978 

1 

1.000 

1 

0.974 

1 

1.000 

1 

0.993 

1 

0.917 

1 

0.986 

1 

Subaverage  %  Meeting  90% 

0.953 

0.786 

1 

1.000 

1 

1.000 

1 

1.000 

1 

0.967 

1 

0.782 

0 

0.970 

1 

0.998 

1 

1.000 

1 

1.000 

1 

1.000 

1 

ISubaverage  %  Meeting  90%  j 

0.974 

0.909 

Average 

Total  %  Meeting  90% 

0.962 

0.840 

Simulation  Results  -  Two-Stage  Heuristic  at  99% 


Attachment  14 


5.475 

7.075 

31.205 

33.906 

6.049 

6.703 

9.161 

17.468 

10.987 

17.164 

57.509 

35.970 

20.977 

5.292 


9.744 

22.467 

7.689 

63.440 

10.928 

11.360 

6.916 

70.470 

3.106 

3.576 

28.399 


Avg  OH  Inv 


3138.040 

596.863 

244.951 

228.775 

263.137 

197.549 

626.275 

752.471 

403.627 

533.637 

416.824 

130.049 

337.304 

1273.430 


305.196 

105.245 

393.471 

49.324 

346.667 

120.490 

167.941 

414.863 

237.255 

667.451 

182.745 


$  OH  Inv 


17180.77 

4222.81 

7643.70 

7756.85 
1591.72 
1324.17 
5737.31 

13144.16 

4434.65 

9159.35 

23971.13 

4677.86 
7075.63 
6738.99 


Subtotal 

114659.08 


2973.83 

2364.54 

3025.40 

3129.08 

3788.38 

1368.77 
1161.48 

29235.40 

736.91 

2386.80 

5189.78 


Subtotal 

55360.36 


Total 

170019.44 


#  Orders 


30 

19 

38 

41 

63 
83 
32 

64 
79 
44 

65 

42 


Subtotal 

747 


14 

11 


Subtotal 

250 


Total 

997 


Fillrate 

Meet  90%  ? 

0.999 

1 

0.969 

1 

1.000 

1 

1.000 

1 

0.968 

1 

0.960 

1 

0.991 

1 

0.997 

1 

1.000 

1 

1.000 

1 

1.000 

1 

1.000 

1 

0.967 

1 

1.000 

1 

Subaverage  %  Meeting  90% 

0.989 

1.000 

1.000 

1 

1.000 

1 

1.000 

1 

1.000 

1 

1.000 

1 

0.914 

1 

1.000 

1 

1.000 

1 

1.000 

1 

1.000 

1 

1.000 

1 

Subaverage  %  Meeting  90% 

0.992 

1.000 

Average  Total  %  Meeting  90% 


Attachment  15 


Proofs  of  Lemmas  and  Theorems 
(as  they  apply  to  the  specific  problem  in  this  project) 

•  Lemma  1:  f3(s,  5 )  <  /3(s  +  1,  5 ) 

Proof:  By  contradiction.  Assume  the  opposite,  that  /?(s,5 )  >  /3(s  +  1,5).  It  follows 
that 


l-p(s,S)<l-(3{s  +  l,S), 

H(s,S)  ^  H{s+l,S) 

M(S—s)n  ^  M(S-s-1)ij.’ 

fl+.S)  .  Hjs+l^S} 

M(S-s-l)+m(S-s-l)  ^  M(S—s—l)  ’ 

M(S  -s-l)*H(s,S)<  [M(5  -  s  -  1)  +  m{S  -  s  -  1)]  *  H(s  +  1,  5 ), 
m(s  -s- 1)  e;:;-1  mo>[£(s  -  j  1 2)  -  l(s  ~  j  1 1»  <  [m<s  -  5  -  p + mts  -  a  -  oi  e^-2  mw)(£(s  ->12)-  ns  ~  >  1  p], 
M(S  —  s  -  l)[L(s  +  1  |  2)  -  L{s  +  1  |  1)]  <  e£T2  m(j)[L(S  -j  |  2)  -  L(5  -  j  |  1)], 
Ef=o-2  +  1  |  2)  -  L(s  +  1  |  1)]  <  Efcr2  m(j)[L(5  -  j  \  2)  -  L(5  -  j  \  1)], 

— )-which  contradicts  the  fact  that 

[L(s  +  1  |  2)  -  L(s  +  1  |  1)]  >  [L(S  -  j  |  2)  -  L(5  -  j  |  1)]  for  j  =  0...5  -  s  -  2. 

Hence  what  we  assumed  is  false  and  /3(s,  S)  <  f3(s  +  1,5). 

•  Lemma  2:  0(s,  5)  <  f3(s  +  1,  5  +  1) 

Proof:  By  contradiction.  Assume  the  opposite,  that  (3(s,S)  >  /3(s  + 1,5  +  1).  It  follows 
that 


1  —  P{s,  5)  <  1  —  /3(s  +  1,5  +  1), 

HjsJS}  ff(s+l,S+l) 

M(S-s)v  ^  M(S+1— s— IV 

H(s,S)  <H(s  +  1,5  +  1), 

Sfcr1  m(j)ms  -  J  |  2)  -  L(S-j  I  1)]  <  Efr0-'  my)[L(S  + 1  -  j  |  2)  -  L(S  + 1  -  j  |  1)], 

— >which  contradicts  the  fact  that.  .  . 

[L(S  -  j  |  2)  -  L(S  -  j  |  1)]  >  [L(S  +  1  -  j  |  2)  -  L(S  +  1  -  j  |  1)]  for  j  =  0...S  -  s  -  1. 
Hence  what  we  assumed  is  false  and  fi(s,  5)  <  /3(s  +  1,5  +  1). 
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•  Lemma  3:  c(s,  S )  <  c(s  +  1,  S ) 

Already  proven  as  Lemma  0  in  [14] . 

•  Lemma  4:  c(s,  S)  <  c(s  +  1,  S  +  1) 

Proof:  By  contradiction.  Assume  the  opposite,  that  c(s,  S)  >  c(s  +  1,  S  +  1).  It  follows 
that 


0  1  m{j)[S-j-2ti+L(S-j\2)]  ^  K+h^I^1  m(j)[S+l-j-2^+L(S+l-j\2)) 

M(S-s)  >  Af(S-s)  > 

Ejsr1  "»(j')[S  -  -  2/1  +  US  -  j  I  2)]  >  ’U,){S  +  1  -  j  -  2m  +  L(S  +  1  -  j  |  2)], 

Efir1  mO')H  +  £(S  -  j  i  2)]  >  Efi’-1  mO)[l  -  i  +  US  +  1  -  j  |  2)], 

Efcr1  m(j)H  +  L(S  -  j  I  2)]  >  Efrr1  mO)[l  -  J  +  L{S  -  j  \  2)  -  F(S  -  j  |  2)], 

where  F(S  -  j  \  2)  =  1  -  E,5  p(j)  <  1. 

— >-The  final  statement  then  contradicts  the  fact  that.  .  . 

H  +  L{S  -  j  I  2)]  <  [1  -  j  +  L(S  -  j  I  2)  -  F(S  -  j  I  2)]  for  j  =  0...S  -s-  1. 

Hence  what  we  assumed  is  false  and  c(s,  S)  <  c(s  +  1,  S'  +  1). 

•  Theorem  1:  If  ( s *,  s*  + 1)  is  the  minimum  feasible  policy  such  that  /3(s*,  s*  + 1)  >  j3, 
then  all  policies  with  S  <  s*  are  infeasible. 

•  Theorem  2:  For  every  S  >  s*  +  1  there  is  one  s  for  which  (a)  all  policies  with  s  <  s 
and  S  <  S  are  infeasible  and  (b)  all  policies  with  s  >  s  and  S  >  S  are  dominated. 

Proofs:  The  proofs  are  by  simple  manipulation  of  the  lemmas  above. 
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Simulation  Results  -  Almost-Exact  Algorithm  at  90%  Attachment  1 7 


c 

Avg  OH  Inv 

$  OH  Inv 

#  Orders 

Fillrate 

Meet  90%? 

A01 

5.475 

1177.255 

6445.47 

78 

0.866 

0 

H mm 

7.075 

235.196 

1664.01 

53 

0.743 

0 

BBS® 

31.205 

114.559 

3574.81 

40 

0.997 

1 

A04 

33.906 

110.294 

3739.63 

21 

0.972 

1 

AOS 

6.049 

94.118 

569.32 

59 

0.789 

0 

A06 

6.703 

69.314 

464.61 

49 

0.803 

0 

A07 

9.161 

274.608 

2515.68 

90 

0.895 

0 

A08 

17.468 

344.431 

6016.52 

73 

0.963 

1 

A09 

10.987 

193.628 

2127.39 

32 

1.000 

1 

A10 

17.164 

262.480 

4505.21 

84 

0.943 

1 

All 

57.509 

214.118 

12313.71 

70 

1.000 

1 

A12 

35.970 

54.284 

1952.60 

44 

0.953 

1 

A13 

20.977 

135.882 

2850.40 

65 

0.851 

0 

A14 

5.292 

539.020 

2852.49 

41 

0.961 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  I 

51591.86 

799 

0.910 

0.571 

B01 

9.744 

138.922 

1353.66 

0.978 

1 

B02 

22.467 

49.755 

1117.85 

1.000 

1 

B03 

7.689 

181.255 

1393.67 

21 

1.000 

1 

B04 

63.440 

20.265 

1285.61 

11 

1.000 

1 

B05 

10.928 

150.000 

1639.20 

31 

0.922 

1 

B06 

11.360 

46.078 

523.45 

43 

0.672 

0 

B07 

6.916 

75.098 

519.38 

13 

0.930 

1 

B08 

70.470 

197.549 

13921.28 

75 

0.982 

1 

B09 

3.106 

115.294 

358.10 

0.944 

1 

BIO 

3.576 

308.039 

1101.55 

0.973 

1 

B12 

28.399 

82.255 

2335.96 

0.967 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  1 

25549.70 

260 

0.943 

0.909 

Total 

Total 

Average 

Total  %  Meeting  90% 

77141.55 

1059 

0.924 

0.720 

Simulation  Results  -  Almost-Exact  Algorithm  at  95%  Attachment  1 7 


c 

Avg  OH  Inv 

$  OH  Inv 

#  Orders 

Fillrate 

Meet  90%? 

A01 

5.475 

1747.549 

9567.83 

78 

0.946 

1 

A02 

7.075 

338.628 

2395.79 

51 

0.852 

0 

A03 

31.205 

154.510 

4821.48 

40 

1.000 

1 

A04 

33.906 

144.510 

4899.76 

21 

1.000 

1 

A05 

6.049 

136.078 

823.14 

61 

0.874 

0 

A06 

6.703 

101.078 

677.53 

49 

0.881 

0 

A07 

9.161 

378.726 

3469.51 

83 

0.941 

1 

A08 

17.468 

470.804 

8224.00 

73 

0.978 

1 

A09 

10.987 

251.961 

2768.30 

40 

1.000 

1 

A10 

17.164 

347.314 

5961.30 

90 

0.971 

1 

All 

57.509 

275.235 

15828.49 

49 

1.000 

1 

A12 

35.970 

75.696 

2722.79 

49 

0.987 

1 

A13 

20.977 

195.177 

4094.23 

72 

0.913 

1 

A14 

5.292 

756.275 

4002.21 

42 

0.986 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  j 

70256.34 

798 

0.952 

0.786 

9.744 

188.333 

1835.12 

27 

1.000 

1 

B02 

22.467 

65.245 

1465.86 

8 

1.000 

1 

7.689 

245.177 

1885.17 

23 

1.000 

1 

B04 

63.440 

28.265 

1793.13 

11 

1.000 

1 

B05 

10.928 

204.314 

2232.74 

36 

0.952 

1 

B06 

11.360 

61.078 

693.85 

46 

0.760 

0 

B07 

6.916 

100.098 

692.28 

16 

0.980 

1 

B08 

70.470 

264.451 

18635.86 

73 

0.997 

1 

B09 

3.106 

154.510 

479.91 

4 

1.000 

1 

BIO 

3.576 

407.451 

1457.04 

14 

1.000 

1 

B12 

28.399 

107.451 

3051.50 

10 

1.000 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  j 

34222.46 

268 

0.972 

0.909 

Total 

Total 

Average 

Total  %  Meeting  90% 

104478.80 

1066 

0.961 

0.840 

Simulation  Results  -  Almost-Exact  Algorithm  at  99%  Attachment  17 

AvgOHInv  $  OH  Inv  |#  Orders  |  Fillrate  Meet90%?  | 


|A01 

5.475 

3139.412 

17188.28 

85 

0.997 

1 

m 

7.075 

575.098 

4068.82 

58 

0.971 

1 

31.205 

244.951 

7643.70 

42 

1.000 

1 

A04 

33.906 

224.510 

7612.24 

21 

1.000 

1 

A05 

6.049 

251.471 

1521.15 

55 

0.973 

1 

A06 

6.703 

174.118 

1167.11 

60 

0.955 

1 

A07 

9.161 

616.078 

5643.89 

81 

0.995 

1 

A08 

17.468 

745.412 

13020.86 

95 

0.996 

1 

A09 

10.987 

396.177 

4352.80 

35 

1.000 

1 

A10 

17.164 

534.755 

9178.53 

90 

1.000 

1 

All 

57.509 

416.118 

23930.53 

70 

1.000 

1 

A12 

35.970 

130.049 

4677.86 

44 

1.000 

1 

A13 

20.977 

334.696 

7020.92 

72 

0.965 

1 

A14 

5.292 

1267.843 

6709.43 

43 

1.000 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90% 

113736.11 

851 

0.989 

1.000 

9.744  305.098  2972.87  34  1.000  1 

22.467  101.520  2280.85  14  1.000  1 

7.689  385.177  2961.63  23  1.000  1 

63.440  48.069  3049.50  13  1.000  1 

10.928  339.118  3705.88  31  1.000  1 

11.360  110.980  1260.73  46  0.910  1 

6.916  164.020  1134.36  13  1.000  1 

70.470  414.157  29185.64  72  1.000  1 

3.106  234.510  728.39  4  1.000  1 

3.576  667.451  2386.80  14  1.000  1 

28.399  182.451  5181.43  10  1.000  1 


Subtotal 

54848.09 

Subtotal 

274 

Subaverage  %  Meeting  90% 

0.992  1.000 

Total 

168584.19 

Total 

1125 

Average  Total  %  Meeting  90% 

0.990  1.000 

Ex  Post  Set  GOF  Test  Results  -  A  Items 


Attachment  18 


A01 


Interval 

Prob 

0 

E 

Stat 

Othr  13 

0.1620 

16 

16.5274 

0.0168 

14  thr  30 

0.1656 

9 

16.8946 

3.6890 

31  thr  52 

0.1680 

13 

17.1365 

0.9985 

53  thr  82 

0.1643 

16 

16.7598 

0.0344 

83  thr  133 

0.1666 

17 

16.9884 

0.0000 

134  thr  inf 

0.1735 

31 

17.6933 

10.0076 

14.7464 

A02 

Interval 

Prob 

0 

E 

Stat 

0  thr  1 

0.1884 

28 

19.2204 

4.0104 

2  thr  3 

0.1300 

3 

13.2644 

7.9429 

4  thr  7 

0.1878 

17 

19.1603 

0.2436 

8  thr  12 

0.1567 

8 

15.9797 

3.9848 

13  thr  21 

0.1630 

15 

16.6283 

0.1594 

22  thr  inf 

0.1740 

31 

17.7470 

9.8970 

26.2381 

A03 

Interval 

Prob 

0 

E 

Stat 

0 

0.2504 

41 

25.5359 

9.3647 

1 

0.1245 

9 

12.6982 

1.0771 

2  thr  3 

0.1576 

19 

16.0720 

0.5334 

4  thr  6 

0.1442 

14 

14.7076 

0.0340 

7  thr  12 

0.1542 

9 

15.7279 

2.8780 

13  thr  inf 

0.1692 

10 

17.2583 

3.0526 

16.9399 

A04 

Interval 

Prob 

0 

E 

Stat 

0 

0.2170 

58 

22.1317 

58.1306 

1 

0.1322 

12 

13.4872 

0.1640 

2  thr  3 

0.1788 

14 

18.2326 

0.9826 

4  thr  6 

0.1662 

5 

16.9491 

8.4241 

7  thr  11 

0.1496 

4 

15.2583 

8.3069 

12  thr  inf 

0.1563 

9 

15.9410 

3.0222 

79.0305 

AOS 

Interval 

Prob 

0 

E 

Stat 

0 

0.2293 

18 

23.3911 

1.2425 

1 

0.1446 

6 

14.7540 

5.1940 

2 

0.1084 

9 

11.0585 

0.3832 

3  thr  4 

0.1548 

13 

15.7853 

0.4915 

5  thr  9 

0.2038 

22 

20.7895 

0.0705 

10  thr  inf 

0.1590 

34 

16.2217 

19.4844 

26.8661 

A06 

Interval 

Prob 

0 

E 

Stat 

Othr  1 

0.2307 

20 

23.5353 

0.5311 

2 

0.1179 

4 

12.0308 

5.3607 

3  thr  4 

0.1993 

7 

20.3296 

8.7399 

5  thr  6 

0.1465 

19 

14.9409 

1.1027 

7  thr  9 

0.1406 

19 

14.3441 

1.5113 

10  thr  inf 

0.1649 

33 

16.8193 

15.5664 

32.8120 

A07 

Interval 

Prob 

0 

E 

Stat 

Othr  6 

0.1645 

15 

16.7808 

0.1890 

7  thr  11 

0.1729 

17 

17.6337 

0.0228 

12  thr  16 

0.1607 

20 

16.3890 

0.7956 

17  thr  23 

0.1780 

23 

18.1600 

1.2900 

24  thr  32 

0.1493 

12 

15.2253 

0.6833 

33  thr  inf 

0.1746 

15 

17.8112 

0.4437 

3.4243 


A08 


Interval 

Prob 

0 

E 

Stat 

Othr  3 

0.1520 

20 

15.5041 

1.3037 

4  thr  8 

0.1980 

13 

20.1971 

2.5646 

9  thr  13 

0.1640 

18 

16.7261 

0.0970 

14  thr  19 

0.1490 

18 

15.2021 

0.5150 

20  thr  30 

0.1695 

22 

17.2869 

1.2850 

31  thr  inf 

0.1675 

11 

17.0838 

2.1665 

7.9318 

A09 

Interval 

Prob 

0 

E 

Stat 

0 

0.1130 

13 

11.5224 

0.1895 

1  thr  2 

0.1998 

27 

20.3770 

2.1527 

3  thr  5 

0.2245 

32 

22.8992 

3.6169 

6  thr  8 

0.1532 

16 

15.6270 

0.0089 

9  thr  13 

0.1523 

9  . 

15.5353 

2.7492 

14  thr  inf 

0.1572 

5 

16.0392 

7.5978 

16.3150 

A10 

Interval 

Prob 

0 

E 

Stat 

Othr  30 

0.1623 

26 

16.5500 

5.3959 

31  thr  59 

0.1642 

23 

16.7500 

2.3321 

60  thr  95 

0.1726 

19 

17.6036 

0.1108 

96  thr  138 

0.1561 

13 

15.9229 

0.5365 

139  thr  21 5 

0.1733 

12 

17.6790 

1.8242 

216  thr  inf 

0.1715 

9 

17.4945 

4.1245 

14.3241 

All 

Interval 

Prob 

0 

E 

Stat 

Othr  5 

0.1720 

40 

17.5399 

28.7604 

6  thr  12 

0.1503 

18 

15.3317 

0.4644 

13  thr  23 

0.1775 

23 

18.1098 

1.3205 

24  thr  38 

0.1665 

8 

16.9828 

4.7514 

39  thr  63 

0.1618 

13 

16.5035 

0.7437 

64  thr  inf 

0.1719 

0 

17.5323 

17.5323 

53.5726 

A12 

Interval 

Prob 

0 

E 

Stat 

0 

0.3710 

30 

37.8426 

1.6253 

1 

0.1760 

18 

17.9502 

0.0001 

2 

0.1139 

12 

11.6150 

0.0128 

3  thr  4 

0.1391 

12 

14.1872 

0.3372 

5  thr  7 

0.1038 

21 

10.5874 

10.2408 

8  thr  inf 

0.0963 

9 

9.8177 

0.0681 

12.2843 

A13 

Interval 

Prob 

0 

E 

Stat 

0  thr  1 

0.1440 

15 

14.6893 

0.0066 

2  thr  4 

0.1901 

15 

19.3927 

0.9950 

5  thr  8 

0.1946 

16 

19.8481 

0.7461 

9  thr  12 

0.1396 

9 

14.2435 

1.9303 

13  thr  20 

0.1689 

26 

17.2230 

4.4729 

21  thr  inf 

0.1628 

21 

16.6034 

1.1642 

9.3150 

A14 

Interval 

Prob 

0 

E 

Stat 

0 

0.1166 

43 

11.8928 

81.3656 

1  thr  4 

0.2172 

7 

22.1505 

10.3627 

5  thr  10 

0.1829 

7 

18.6559 

7.2824 

11  thr  17 

0.1358 

6 

13.8537 

4.4522 

18  thr  34 

0.1795 

10 

18.3070 

3.7694 

35  thr  inf 

0.1680 

29 

17.1402 

8.2061 

115.4383 


Ex  Post  Set  GOF  Test  Results  -  B  Items 


Attachment  18 


B01 

Interval 

Prob 

o 

E 

Stat 

B07 

Interval 

Prob 

O 

E 

Stat 

0 

0.2634 

35 

26.8650 

2.4634 

0 

0.5737 

65 

58.5201 

0.7175 

1 

0.1564 

21 

15.9487 

1.5999 

1 

0.1626 

14 

16.5852 

0.4030 

2 

0.1131 

13 

11.5356 

0.1859 

2 

0.0869 

14 

8.8655 

2.9737 

3  thr  4 

0.1550 

12 

15.8105 

0.9184 

3 

0.0537 

2 

5.4811 

2.2109 

5  thr  8 

0.1649 

14 

16.8159 

0.4715 

4  thr  5 

0.0598 

3 

6.0973 

1.5733 

9  thr  inf 

0.1473 

7 

15.0243 

4.2857 

6  thr  inf 

0.0632 

4 

6.4509 

0.9312 

9.9247 

8.8096 

B02 

B08 

Interval 

Prob 

0 

E 

Stat 

Interval 

Prob 

0 

E 

Stat 

0 

0.4584 

75 

46.7530 

17.0661 

27 

16.6564 

6.4234 

1 

0.1780 

13 

18.1561 

1.4643 

6  thr  14 

0.1701 

15 

17.3477 

0.3177 

2 

0.1057 

6 

10.7805 

2.1199 

15  thr  26 

0.1654 

20 

16.8758 

0.5784 

3 

0.0700 

5 

7.1393 

0.6410 

27  thr  44 

0.1693 

14 

17.2722 

0.6199 

4  thr  5 

0.0837 

2 

8.5377 

5.0062 

45  thr  74 

0.1617 

18 

16.4962 

0.1371 

6  thr  inf 

0.1042 

1 

10.6334 

8.7274 

75  thr  inf 

0.1701 

8 

17.3517 

5.0401 

35.0248 

13.1166 

B03 

B09 

Interval 

Prob 

0 

E 

Stat 

Interval 

Prob 

0 

E 

Stat 

0 

49 

35.7826 

0 

0.8623 

91 

87.9525 

1 

0.1493 

18 

15.2304 

0.5036 

1 

0.0952 

6 

9.7088 

1.4168 

2 

0.0980 

8 

9.9991 

0.3997 

2 

0.0273 

4 

2.7810 

0.5343 

3  thr  4 

0.1274 

9 

12.9958 

1.2286 

3 

0.0094 

0 

0.9598 

0.9598 

5  thr  8 

0.1336 

13 

13.6322 

0.0293 

4 

0.0035 

1 

0.3594 

1.1417 

9  thr  inf 

0.1408 

5 

14.3600 

6.1009 

5  thr inf 

0.0023 

0 

0.2384 

0.2384 

13.1444 

4.3965 

B04 

BIO 

Interval 

Prob 

0 

E 

Stat 

Interval 

Prob 

0 

E 

Stat 

0 

0.7573 

79 

77.2441 

0.0399 

0 

0.6229 

82 

63.5351 

5.3664 

1 

0.0739 

3 

2.7337 

1 

0.1189 

0 

12.1267 

12.1267 

2 

0.0379 

0 

3.8706 

3.8706 

2 

0.0643 

6 

6.5537 

0.0468 

3  thr  4 

0.0426 

12 

4.3431 

13.4990 

3  thr  4 

0.0724 

4 

7.3850 

1.5516 

5  thr  9 

0.0463 

3 

4.7219 

0.6279 

5  thr  7 

0.0540 

5 

5.5089 

0.0470 

10  thr  inf 

0.0420 

5 

4.2801 

0.1211 

8  thr  inf 

0.0676 

5 

6.8906 

0.5187 

20.8923 

M 

19.6572 

B05 

B12 

Interval 

Prob 

O 

E 

Stat 

Interval 

Prob 

O 

E  ' 

Stat 

0 

36 

29.4604 

1.4516 

0 

0.8005 

89 

81.6532 

1 

0.1493 

8 

15.2254 

3.4289 

1 

0.0680 

2 

6.9387 

3.5152 

2 

0.1042 

9 

10.6263 

0.2489 

2 

0.0337 

0 

3.4414 

3.4414 

3  thr  4 

0.1413 

12 

14.4117 

0.4036 

3  thr  4 

0.0363 

2 

3.7062 

0.7854 

5  thr  8 

0.1538 

20 

15.6858 

1.1865 

5  thr  8 

0.0319 

3 

3.2525 

0.0196 

9  thr  inf 

0.1627 

17 

16.5904 

0.0101 

9  thr  inf 

0.0295 

6 

3.0080 

2.9760 

6.7296 

11.3987 

B06 

Interval 

Prob 

0 

E 

Stat 

0 

0.3544 

39 

36.1501 

0.2247 

1 

0.1929 

12 

19.6786 

2.9962 

2 

0.1267 

7 

12.9226 

2.7144 

3 

0.0879 

3 

8.9699 

3.9733 

4  thr  5 

0.1081 

12 

11.0240 

0.0864 

6  thr  inf 

0.1299 

29 

13.2549 

18.7033 

28.6982 


Simulation  Results  -  Almost-Exact/Croston  Hybrid  Attachment  19 


c 

Avg  OH  Inv 

$  OH  Inv 

#  Orders 

Fillrate 

Meet  90%? 

A01 

5.475 

1177.255 

6445.47 

78 

0.866 

0 

mm- 

7.075 

549.314 

3886.40 

54 

0.957 

1 

BS& 

31.205 

73.480 

2292.94 

43 

0.930 

1 

A04 

33.906 

102.206 

3465.40 

23 

0.975 

1 

A05 

6.049 

257.843 

1559.69 

62 

0.994 

1 

A06 

6.703 

69.314 

464.61 

49 

0.803 

0 

A07 

9.161 

274.608 

2515.68 

90 

,0.895 

0 

A08 

17.468 

344.431 

6016.52 

73 

0.963 

1 

A09 

10.987 

193.628 

2127.39 

32 

1.000 

1 

A10 

17.164 

262.480 

4505.21 

84 

0.943 

1 

All 

57.509 

104.529 

6011.36 

54 

0.974 

1 

A12 

35.970 

54.284 

1952.60 

44 

0.953 

1 

A13 

20.977 

135.882 

2850.40 

65 

0.851 

0 

A14 

5.292 

762.255 

4033.85 

46 

0.984 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90% 

48127.52 

797 

0.935 

0.714 

B01 

9.744 

138.922 

1353.66 

27 

1 

B02 

22.467 

20.539 

461.45 

20 

0.963 

1 

B03 

7.689 

181.255 

1393.67 

21 

1.000 

1 

B04 

63.440 

‘  5.21-6 

330.90 

21 

0.705 

0 

B05 

10.928 

119.706 

1308.15 

52 

0 

B06 

11.360 

82.157 

933.30 

49 

0.833 

0 

B07 

6.916 

75.098 

519.38 

13 

0.930 

1 

B08 

70.470 

148.235 

10446.12 

67 

0.941 

1 

B09 

3.106 

115.294 

358.10 

4 

0.944 

1 

BIO 

3.576 

108.628 

388.45 

20 

0.761 

0 

B12 

28.399 

82.255 

2335.96 

9 

0.967 

1 

Subtotal 

Subtotal 

Subaverage  %  Meeting  90%  J 

19829.14 

303 

0.901 

0.636 

Total 

Total 

Average 

Total  %  Meeting  90% 

67956.66 

1100 

0.920 

0.680 

Designing  an  (s,  5)  Inventory  Contror  System 

Gary  Rafnson 
18  July  1997 

Abstract  :  ' is 'i'  '  ‘ 

v  This  project  aims  to  provide  Borne  insight  on  a  local  company’s  ($,  S)  inventory 
control  system.  It  does  so  by  providing  three  completely  |£W  (a,  S)  approaches  and 
comparing  their  simulated  performance  to  the  current  system]'  Recommendations  are 
made  based  on  the  results.  !  V 
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Abstract 
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